Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,759,226 | 1 | . | G | 66.7% | 45.2 / 26.6 | 24 | intergenic (‑91/‑51) | yfjM/rnlA | DUF1524 domain‑containing protein/type II toxin‑antitoxin system toxin mRNA endoribonuclease RnlA |
Reads supporting (aligned to +/- strand): ref base . (7/1); new base G (8/8); total (15/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCA > NZ_CP009273/2759137‑2759313 | ttGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCa > 2:234868/1‑90 (MQ=255) gtgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑gg < 1:240943/90‑1 (MQ=255) gACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑ggg‑gg < 2:762917/90‑1 (MQ=255) tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑ggg‑ggggt < 2:414865/90‑3 (MQ=255) gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑ggg‑ggggta > 2:733645/1‑87 (MQ=255) cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑ggg‑ggggtat < 2:775066/90‑5 (MQ=255) cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑ggg‑ggggtat < 1:821549/90‑5 (MQ=255) cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑GGGCGGGtat < 1:642123/90‑1 (MQ=255) cccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTattatt < 2:315861/90‑1 (MQ=255) aCGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGtt < 1:551183/90‑1 (MQ=255) aCGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGtt < 1:234868/90‑1 (MQ=255) aCGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGtt < 2:21110/90‑1 (MQ=255) gCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTCCCCata > 2:258886/1‑90 (MQ=255) gCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCatat > 2:151572/1‑90 (MQ=255) ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTCCCCATATGttt > 1:830912/1‑90 (MQ=255) ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTCCCCATATGttt > 2:491069/1‑90 (MQ=255) tGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATgg > 2:796280/1‑90 (MQ=255) cggAGCGCCCCCTATTCAGGCGCCTAAATTAAGCAAATCACTACACTTCAGGGG‑GGGTATTATTGTAGTGTTTCCCCATATGTTTCTATg > 2:172935/2‑90 (MQ=255) aGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTccc > 1:588940/1‑73 (MQ=255) aGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTccc < 2:588940/73‑1 (MQ=255) gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCa > 2:476504/1‑90 (MQ=255) gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCa > 2:791428/1‑90 (MQ=255) cAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACa > 2:190789/1‑90 (MQ=255) ggTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTccc > 2:348061/1‑58 (MQ=255) ggTGGCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTccc < 1:348061/58‑1 (MQ=255) ggCTAAATTTAGCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa > 1:320219/1‑90 (MQ=255) gCAAATCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaa < 2:79500/90‑1 (MQ=255) aaTCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGa > 2:484557/1‑58 (MQ=255) aaTCAATACACTTCAGGGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGa < 1:484557/58‑1 (MQ=255) cTTCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGt > 2:560212/1‑90 (MQ=255) tCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCa > 2:533653/1‑90 (MQ=255) | TTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCA‑GGG‑GGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCA > NZ_CP009273/2759137‑2759313 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |