| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,436,061 | 0 | T | G | 54.3% | 41.9 / 49.8 | 35 | G273G (GGT→GGG) | mnmC | bifunctional tRNA (5‑methylaminomethyl‑2‑thiouridine)(34)‑methyltransferase MnmD/FAD‑dependent 5‑carboxymethylaminomethyl‑2‑thiouridine(34) oxidoreductase MnmC |
| Reads supporting (aligned to +/- strand): ref base T (3/13); new base G (19/0); total (22/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.79e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGTGC > NZ_CP009273/2435979‑2436150 | aGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGccc > 1:1276321/1‑89 (MQ=255) gACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAg < 2:1537452/90‑1 (MQ=255) tACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCgcgccg > 2:317983/1‑88 (MQ=255) gCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGGgttttt > 1:1062945/1‑90 (MQ=255) cTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCACGCggtgttt > 1:1318774/1‑89 (MQ=255) ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCCGCGCGGTGTTGTc > 1:287959/1‑90 (MQ=255) ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCGCGCGGTGTTTTc > 1:1505160/1‑90 (MQ=255) ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGGCAGCGCGGtgttgtg > 2:1401439/1‑89 (MQ=255) ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCCGCGCGGtgttgtt > 1:510904/1‑89 (MQ=255) gCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCCGCGCGGTGTTGTTGCTgg > 1:584339/1‑90 (MQ=255) ccGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCAAGGCGTTGTTGTTGCTGgcgc > 1:1326995/1‑90 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAg > 1:119952/1‑68 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAg < 2:119952/68‑1 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTa > 1:447422/1‑90 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCAGCGCGTTGTTGTCGCTGGCGCTa > 1:1375263/1‑90 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGCGGTGTTGTCGCTGGCGCtt > 1:100422/1‑89 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGGTTGCCCGCGCGGTGGTGTCGCTGgggccg > 2:726231/1‑88 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCGCGCGGTGTTGTTGCTGGCGCTa > 1:1074976/1‑90 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCGCGCGGTGTTGTCGCTGGGGCTa > 2:431459/1‑90 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCCCGCGGTGTTTTCGCCGGCGGTa > 2:641298/1‑90 (MQ=255) gTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCCGCGCGTTGTTGTCGCTGGCGCTATTTCg > 1:925352/1‑90 (MQ=255) gTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGGTTTCCCGGGCGGTGTTGTTGCTGGGGCTAtttag > 2:1223913/1‑87 (MQ=255) ccGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGTTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCggagc < 1:1085432/90‑4 (MQ=255) cACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTgggg > 1:226099/1‑88 (MQ=255) aCGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCCGCGCGGTGTTGTCGCTGGCGCTATTACGGCGCGGCTgggcg > 1:763078/1‑87 (MQ=255) gggCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGt < 1:933739/90‑1 (MQ=255) gggCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGt < 2:925352/90‑1 (MQ=255) gcagcagcaAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTaa < 1:431459/90‑1 (MQ=255) aaaCGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTAt < 1:130670/90‑1 (MQ=255) ggAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTgcgc < 1:1223913/90‑1 (MQ=255) ggAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTgcgc < 2:763078/90‑1 (MQ=255) gcgATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAg < 2:100422/90‑1 (MQ=255) gcgATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAg < 2:447422/90‑1 (MQ=255) aTCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGccccc < 2:1326995/90‑1 (MQ=255) aTCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGccccc < 2:499340/90‑1 (MQ=255) gtggtATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTgg < 2:287959/90‑1 (MQ=255) tggtATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTggg < 1:1589782/90‑1 (MQ=255) tATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGTGc < 2:1276321/90‑1 (MQ=255) | AGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGGTGC > NZ_CP009273/2435979‑2436150 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |