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% | 25.5 / 37.8 | 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.62e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC > NZ_CP009273/2778606‑2778774 | aaTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc < 1:159083/90‑1 (MQ=255) aaTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc < 2:194395/90‑1 (MQ=255) aTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa < 1:376483/90‑1 (MQ=255) ctcAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGa < 1:249448/90‑1 (MQ=255) cAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTTGTCGGTAGTAGATg > 1:824937/1‑90 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc > 2:829332/1‑72 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc < 1:829332/72‑1 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCAGTAGGAGGTGttttaag > 2:832335/1‑86 (MQ=255) ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAGATGTTTTAAGCGt > 2:392234/1‑90 (MQ=255) gAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCagag < 1:698449/90‑1 (MQ=255) gAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCagag < 1:832335/90‑1 (MQ=255) aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGTCCGTAGTAGAAGTTTAAGGCGGGGGAGAGACATTTCATc > 1:468942/1‑90 (MQ=255) gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTTGTCAGTAGGAGATGTTTTAAGCGGGGCAGAGAGATTTTATCCt > 1:493702/1‑90 (MQ=255) aaCGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTAc < 2:468942/90‑1 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc < 1:587379/90‑1 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAGATGTTTTAAGCGTGGGAGAGACATTTTATCCCTTCTCTACGGc > 1:246736/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCCGGAGTAGATGTTTAAAGGGTGGGCGAGACAATTTATTCTTTCTCTACCGc > 1:383487/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTTGGCAGTAGTAGATGTTTAAGGCGGGGGAGAGAAATTTTTTCCTTTCTCTACGGc > 2:828714/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc > 2:163708/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTTATCCTTTCTCTACGGc > 2:436692/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGGAGATGTTTAAGGCGGGGCAGAGACATTTTATCCTTACTCTACGGc > 2:623821/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGAAGTTTAAAGGGGGGGGGGGGAAATTCATCCTTTCTCTACCGc > 1:429314/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGTCCGTAGTAGATGTTTAAGGGGGGGCAGAGACATTTTATCCTTTCTCTTCGGc > 1:781059/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 2:641919/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTTATCCTTTCTCTACGGCAt > 1:316533/1‑90 (MQ=255) cacaAAACTCATCTGCCTGCGGCGGGGGTTTTTGTCAGTAGTAGATGTTTTAGGCGGGGGAGGGACATTTTATCCTTTCTCTACGGCAtt > 1:88898/1‑90 (MQ=255) cATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacat < 1:623821/90‑1 (MQ=255) gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc < 2:412334/90‑1 (MQ=255) gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc < 2:429314/90‑1 (MQ=255) | AATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC > NZ_CP009273/2778606‑2778774 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |