Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,778,686 | 0 | T | G | 55.6% | ‑3.6 / 27.9 | 27 | intergenic (‑316/+435) | BW25113_RS25705/BW25113_RS13825 | DUF5507 domain‑containing protein/tRNA‑Ile |
Reads supporting (aligned to +/- strand): ref base T (1/11); new base G (15/0); total (16/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.20e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC > NZ_CP009273/2778601‑2778774 | ttAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTc < 2:777304/90‑1 (MQ=255) ctcAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGa < 2:279772/90‑1 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc < 2:602128/72‑1 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa < 2:306460/73‑1 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTgggc > 1:602128/1‑70 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCa > 1:306460/1‑73 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGGAGGAAAAGttttaag > 1:530692/1‑86 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTAAgg > 2:201424/1‑90 (MQ=255) tGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCg < 2:157436/90‑1 (MQ=255) ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCTTCTGCCTGCGGCGGGGGTTCTGGACAGAAGGAGATGGTTTAGGGGt > 1:341449/1‑90 (MQ=255) ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc < 2:249630/69‑1 (MQ=255) ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTgggc > 1:249630/1‑67 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGGAGAGACa > 1:329904/1‑90 (MQ=255) gCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTc < 2:652196/90‑1 (MQ=255) gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTTAAGCGTGGCCGAGACATTTTATCCt > 1:381569/1‑90 (MQ=255) gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTTGTCCGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTTATCCt > 2:459240/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGGAGATGTTTTAGGCGGGGGCGAGACAATTTATCCTTTCTCTTCGGc > 2:251824/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCCCTAATAGATGTTTAAGGCGGGGGAAAGACATTTCATTCTTACTCTACGGc > 1:363506/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTTATCCTTACTCTACGGc > 1:690739/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 2:238802/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 2:329904/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 2:363506/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTTATCCTTTCTCTACCGCAt > 1:642383/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGGAGTAGATGTTTTAGGCGGGGCAGAGACATTTTATCCTTTCTCTTCGGCAt > 1:434695/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGTCAGTAGTAGATGTTTTAGGCGTGGCAGAGACATTTTATCCTTTCTCTACGGCAt > 2:363822/1‑90 (MQ=255) ctgcctgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacataca < 1:579610/90‑1 (MQ=255) gcctgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACAtt > 1:128354/1‑90 (MQ=255) gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc < 2:100073/90‑1 (MQ=255) | TTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC > NZ_CP009273/2778601‑2778774 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |