Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,450,688 | 0 | T | G | 53.3% | 10.4 / 17.6 | 15 | Y597S (TAC→TCC) | nrdD | anaerobic ribonucleoside‑triphosphate reductase |
Reads supporting (aligned to +/- strand): ref base T (0/7); new base G (8/0); total (8/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCCAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGC‑GCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGAAACTGTTGGTGTAGTAACCTTT > NZ_CP009273/4450601‑4450770 | tCCCAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTaa < 2:4745/90‑1 (MQ=255) cccAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGCTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGAAg > 1:57946/1‑90 (MQ=255) cAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGgc < 1:323412/90‑1 (MQ=255) aGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGAAgggg > 2:215252/1‑88 (MQ=255) cAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGC‑GCTTCAAAGTc < 1:74849/90‑1 (MQ=255) ttCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGAAGGC‑GGGTCAAAGACGGTCttt > 1:550639/1‑89 (MQ=255) gtgCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGC‑GCTTCAAAGTCGATCTTGTCGTACgg < 1:69396/90‑1 (MQ=255) gtgCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGGAAGG‑GCCTCCAAATCGATCTTTTCGTTCgg > 2:225569/1‑90 (MQ=255) gtgCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGAAGGG‑GGCTCCAAAACGATTTTGTCGTACgg > 2:96869/1‑90 (MQ=255) tttGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGAAGGG‑GCTTCAAAATCGATCTTTTCGTTCGGGGTTACcttc > 2:69396/1‑90 (MQ=255) cTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGGAGGCGGGTTGAAAAT‑GATCTTGTCGTACGGGGTTCCCCTCTTTTcccc > 2:205296/1‑88 (MQ=255) aaTGAAACCACCGTTCGCCAGCGGCGGGTAAGGC‑GCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGaa < 1:545095/90‑1 (MQ=255) accaccGTTCGCCAGCGGCGGGTAAGGC‑GCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGAAACTGtt < 1:216547/90‑1 (MQ=255) ccGTTCGCCAGCGGCGGGGAAAGG‑GCCTCAAAATCGATCTTTTTGTACGGGGTCACCCTCTTCTCCCCATCCAGGGGGGAAATGTTGgtg > 1:490313/1‑90 (MQ=255) cggcggGTAAGGC‑GCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGAAACTGTTGGTGTAGTAACCttt < 1:518703/90‑1 (MQ=255) | TCCCAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGC‑GCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGAAACTGTTGGTGTAGTAACCTTT > NZ_CP009273/4450601‑4450770 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |