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.8% | 9.1 / 15.4 | 13 | Y597S (TAC→TCC) | nrdD | anaerobic ribonucleoside‑triphosphate reductase |
Reads supporting (aligned to +/- strand): ref base T (0/6); new base G (7/0); total (7/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCCAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGCGCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGAAACTGTTGGTGTAGTAA > NZ_CP009273/4450601‑4450765 | tCCCAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTaa < 1:220025/90‑1 (MQ=255) cAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGAAggg > 1:312370/1‑89 (MQ=255) ttGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGCGCTTCAAAGTCGATCTTGTCGTAc < 2:212644/90‑1 (MQ=255) gtgCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGGAGGGGCCTCCAAAACGATCTTGTCGTACgg > 2:168962/1‑90 (MQ=255) gtgCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGAAGGGGCTTCAAAAACGATCTTTTCGTACgg > 1:356458/1‑90 (MQ=255) tttGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGCGCTTCAAAGTCGATCTTGTCGTACGGGTTCACcttc < 1:42041/90‑1 (MQ=255) aCTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGGAGGGGCGTCTAAAACGAATTTTTCGTTCGGGGTTACCCTCTTTTccc > 2:341898/1‑89 (MQ=255) gCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGCGCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATc < 2:267027/90‑1 (MQ=255) gTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGGAAGGCGCTTCAAAGGCGATCTTGTCGTACGGGGTCACCTTCTTCTCCCCATCGAg > 1:50373/1‑90 (MQ=255) ccGTTCGCCAGCGGCGGGGAAAGGGCTTTAAAAACGATCTTGTCGTACGGGTTTACCTTCTTTTCCCCATCGAGGGGGGAAATGTTtgtg > 2:3068/1‑90 (MQ=255) gTTCGCCAGCGGCGGGGAAGGGGCTTCAAAAACGATCTTTTCGTACGGGGTCACCTTCTTTTCCCCATCCAGGGGGAAACTGTTGGGgta > 1:212507/1‑90 (MQ=255) tCGCCAGCGGCGGGTAAGGCGCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGAAACTGTTGGTgtagt < 1:39282/90‑1 (MQ=255) gCCAGCGGCGGGTAAGGCGCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGAAACTGTTGGTGTAGTaa < 1:341898/90‑1 (MQ=255) | TCCCAGACATCTTCCAGCGCCTTCAGGTTGTGCTGAATGTTTGGATACTCGCCGTAGCAAATGAAACCACCGTTCGCCAGCGGCGGGTAAGGCGCTTCAAAGTCGATCTTGTCGTACGGGTTCACCTTCTTCTCCACATCGAGGTGGAAACTGTTGGTGTAGTAA > NZ_CP009273/4450601‑4450765 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |