| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,026,091 | 0 | T | G | 52.6% | 14.7 / 23.6 | 19 | G80G (GGT→GGG) | yigZ | IMPACT family protein |
| Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (10/0); total (11/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCAT > NZ_CP009273/4026004‑4026168 | ccATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTaa < 1:262418/90‑1 (MQ=255) gtggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATg > 1:110779/1‑90 (MQ=255) tggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGc < 1:197495/90‑1 (MQ=255) gcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGGAAACCATGCTc > 2:75156/1‑90 (MQ=255) gtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAAACCATGCTCGc > 2:168817/1‑90 (MQ=255) ggTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGcccc > 1:119501/1‑89 (MQ=255) ggTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCa < 1:376540/90‑1 (MQ=255) ggTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCa < 2:411440/90‑1 (MQ=255) ggTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCa < 1:88038/90‑1 (MQ=255) ggTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCa < 1:168817/90‑1 (MQ=255) gcgGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTa > 1:172817/1‑90 (MQ=255) cAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCCATGCTCGCCCAACTAATGGGCAGCGGCGGCGGGGaa > 1:31274/1‑90 (MQ=255) cAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCCCGCCCCGCCAATGGGGAGCGGGGGCGGGGgaaat > 1:363832/1‑88 (MQ=255) cTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTAcc < 2:87964/90‑1 (MQ=255) gggTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAAACGATGCTCCCCCCGCTAAAGGGGGGCGGGGGCGGGGGAAATAACGc > 2:271702/1‑90 (MQ=255) gacgacGGGGAGCCGGCGGGAACGGCAGGGAAACCCATGCTCGCCCCGCTAATGGGCAGCGGCGTCGGGGGAAATACCGCTGTGGGAGTg > 1:256028/1‑90 (MQ=255) ggggAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCCAATGGGCCGCGGGGGCGGGGGAAATACCGCCGTGGGAGGGCGctac > 2:40428/1‑90 (MQ=255) ggcggGAACGGCAGGGAAAACCATGCTCCCCCCGCTAAAGGGGGGCGGGGGCGGGGGAAATACCCCCGTGGGAGGGGGCCAATAcggcgg > 2:455458/1‑90 (MQ=255) gggAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCat < 1:466877/90‑1 (MQ=255) | CCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCAT > NZ_CP009273/4026004‑4026168 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |