| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,778,686 | 0 | T | G | 42.6% | 76.5 / 10.9 | 69 | intergenic (‑316/+435) | BW25113_RS25705/BW25113_RS13825 | DUF5507 domain‑containing protein/tRNA‑Ile |
| Reads supporting (aligned to +/- strand): ref base T (7/32); new base G (29/0); total (37/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.02e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC > NZ_CP009273/2778599‑2778774 | tgtTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGgtgt < 2:482559/90‑1 (MQ=255) aTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa < 1:92744/90‑1 (MQ=255) ctcAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGa < 2:1187745/90‑1 (MQ=255) cAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTTAGGAGTAAATg > 2:1302954/1‑90 (MQ=255) gCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTaa < 1:1206111/90‑1 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAgg < 2:583162/90‑1 (MQ=255) gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGTGAGTAGTAGATGTTTTAgg > 1:786175/1‑90 (MQ=255) ttGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGGAGGAGATGTTTAAGGc > 2:1206111/1‑90 (MQ=255) ttAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCCGTAGTAGATGTTTTAGGCGTgg > 1:1011275/1‑90 (MQ=255) aTGATCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa < 2:237833/65‑1 (MQ=255) aTGATCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa > 1:237833/1‑65 (MQ=255) aGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAgaga < 1:239324/90‑1 (MQ=255) gCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTGGGAGATGTTTTAAGCGGGGCGGGGAc > 1:200723/1‑90 (MQ=255) gCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAGGCGGGGCAGAGAc > 2:654995/1‑90 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCagtag > 2:358864/1‑64 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCagtag < 1:358864/64‑1 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa > 1:1045738/1‑60 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa < 2:1045738/60‑1 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACa > 2:824175/1‑90 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACa > 2:1235350/1‑90 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGGCAGTAGTAGATGTTTAAGGCGTGGCAGAGACa > 2:798303/1‑90 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAGGCGGGGGAGAGACa > 2:783225/1‑90 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGGGAGTGATAGATGTTTTAATCGTGGCAGAGACa > 2:1110915/1‑90 (MQ=255) cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGGCAGTAGGAGAAGTTTAAAGCGTGGCAGAGACa > 1:1415942/1‑90 (MQ=255) aTAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTTGTTAGTAGGAGATGTTTTAGGGGGGGCGGAGACAt > 1:236118/1‑90 (MQ=255) tAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGACGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACAtt < 1:1011839/90‑1 (MQ=255) gCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTc < 1:1235350/90‑1 (MQ=255) gCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCAGTAGTAGATGTTTTAAGCGGGGGAGAGACAtttt > 1:509432/1‑89 (MQ=255) ccACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACAtttta > 2:831300/1‑88 (MQ=255) aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAATAAGAGATGTTTTAGGCGTGGCGGAGGCCTTTCCTc > 2:1133608/1‑90 (MQ=255) cGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGGCAGTAGTAGATGTTTTAAGGGGGGGGGAGACATTTctttc > 1:1163901/1‑86 (MQ=255) gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTgg < 2:1298502/48‑1 (MQ=255) gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCt < 1:1302954/90‑1 (MQ=255) gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCAGTAGTAGATGTTTAAAGCGTGGCAGAGAAAATTTATCCt > 1:1415797/1‑90 (MQ=255) gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTgg > 1:1298502/1‑48 (MQ=255) cTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCtt < 2:1415942/90‑1 (MQ=255) ccTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTc < 2:179816/42‑1 (MQ=255) ccTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTAc < 1:133631/90‑1 (MQ=255) ccTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTc > 1:179816/1‑42 (MQ=255) aaCGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGGAGAGACATTTCATCCTTACTCTAc < 1:776061/90‑1 (MQ=255) aaCGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTAc < 1:831300/90‑1 (MQ=255) aaCGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTTATCCTTTCTCTAc > 2:523794/1‑90 (MQ=255) aCGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGGAGTAGATGTTTTAAGCGGGGGAGAGACATTTTCTCCTTTCTCTACg > 2:719697/1‑90 (MQ=255) aCGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAATAGTAGATGTTTAAGGGGGGGGGGAGACATTTTATCCCTTCTCTTCg > 2:786103/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc < 1:369891/90‑1 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTCATCCTTACTCTACGGc > 2:502864/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAGATGTTTTAGGGGGGGGAGAGACATTTCATCCTTTCTCTACGGc > 1:694067/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAGGCGGGGGGGAGACATTTTATCCTTTCTCTACGGc > 2:776061/1‑90 (MQ=255) gCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGTCAGTAGGAGATGTTTAAAGCGGGGGAGAGACATTTTATCCTTACTCTTCGGc > 1:639392/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 2:786175/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt < 1:523794/90‑1 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGAAGTTTAAAGCGGGGCAGAGGCATTTTATTCTTTCTCTACGGCAt > 1:22079/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTAAGGCGGGGGAGAGACATTTCATCCTTTCTCTACGGCAt > 2:1416829/1‑90 (MQ=255) tCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGTCAGTAGTAGATGTTTTAGGGGGGGGAGGGAAATTTTATCCTTACTCTACGGCAt > 2:1098238/1‑90 (MQ=255) aaaaCTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGtt < 2:200723/90‑1 (MQ=255) aaaCTCATCTGCCTGCGGCGGGTGTTCTGGTCTGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTc < 2:646155/90‑1 (MQ=255) aaaCTCATCTGCCTGCGGCGGGGGGTTTGGTCCGGAGGAGATGTTTAAGGCGGGGGCGAGACATTTTATCCTTACTCTTCTGCATTGttt > 1:314827/1‑89 (MQ=255) aaCTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCt < 1:714818/90‑1 (MQ=255) cTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtac < 1:502864/90‑1 (MQ=255) tCATCTGCCTGCGGCGGGCGTTCTGGGCAGTAGTAGATGTTTTAGGTGGGGCAGAGACATTTTCTCCTTATTCTACGGCATTGTTTtaca > 2:881413/1‑90 (MQ=255) cATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacat < 1:47907/90‑1 (MQ=255) cATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacat < 1:654995/90‑1 (MQ=255) aTCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacata > 2:1251513/1‑90 (MQ=255) cctgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTg < 1:1251513/90‑1 (MQ=255) tgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGt < 2:1415797/90‑1 (MQ=255) tgcGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGt > 2:525707/1‑90 (MQ=255) cggcggGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTtg < 1:653139/90‑1 (MQ=255) cggcggGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTtg < 2:1163901/90‑1 (MQ=255) cggGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTgg < 1:824175/90‑1 (MQ=255) ggTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTa < 2:972436/90‑1 (MQ=255) gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc < 1:1183501/90‑1 (MQ=255) gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc < 1:881413/90‑1 (MQ=255) | TGTTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC > NZ_CP009273/2778599‑2778774 |
| 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: ‑ |