| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 5,162 | 0 | T | A | 57.1% | 0.7 / 13.4 | 14 | intergenic (+142/‑72) | thrC/yaaX | threonine synthase/DUF2502 domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base T (5/1); new base A (0/8); total (5/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.00e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAA > NZ_CP009273/5079‑5248 | gaagaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATaacaa‑c > 2:282696/1‑90 (MQ=255) gaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCg > 1:168622/1‑90 (MQ=255) aaTGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAATTT‑AGAGATTAGGTTTGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTa < 2:550087/90‑1 (MQ=255) aaGGAAAAATCCAAAAAAAATG‑GGTAAATTAAGAGATTAGGATGGGGAAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATc < 1:247224/90‑1 (MQ=255) aTAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa > 1:157479/1‑90 (MQ=255) aaaaTGGGGTAACTT‑AGAGATTAGGATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 2:157479/89‑1 (MQ=255) aaaaGGCGATAACTT‑AGAGTTTAGGATGGGGAAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 1:292178/89‑1 (MQ=255) ggTAATTA‑AGATATTAGGTTGGGGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 2:337387/90‑1 (MQ=255) ggTAATTA‑AGAAATT‑GGGTTGGGGAAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 1:197881/90‑1 (MQ=255) ggTAACTT‑AGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa > 2:564068/1‑90 (MQ=255) agaTTAGGATGGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc < 2:90193/90‑1 (MQ=255) gaTTAGGATGGCGAAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt < 1:424456/90‑1 (MQ=255) ttAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt > 1:90193/1‑90 (MQ=255) aaaTAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa < 1:202008/89‑1 (MQ=255) | GAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAA > NZ_CP009273/5079‑5248 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |