Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,436,061 | 0 | T | G | 55.6% | 14.6 / 24.0 | 18 | 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 (0/8); new base G (10/0); total (10/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.29e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGG > NZ_CP009273/2435974‑2436147 | ggAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTAt < 1:256054/90‑1 (MQ=255) cAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGcc < 1:582440/90‑1 (MQ=255) aGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGGTTGccc > 2:380636/1‑89 (MQ=255) aGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCa > 1:119723/1‑90 (MQ=255) tACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCcgcgcgg > 1:50818/1‑89 (MQ=255) cTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGCGttggtg > 2:379293/1‑90 (MQ=255) ccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCGCGCGGTGTTGTCGc > 1:73870/1‑90 (MQ=255) ccGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGgcgc < 1:449054/90‑1 (MQ=255) ccGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGgcgc < 1:467261/90‑1 (MQ=255) ccGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTGCCGCGCGGGGTTTTCGCCTgcgc > 1:570062/1‑90 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCCGCGCGGTGTTGTTGCTGgcgcct > 1:133255/1‑88 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGGGTGCCAGCGCGGTGTTGTTGGTGGCGCTa > 2:467509/1‑90 (MQ=255) gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCGTTGTTGTCGCTGGCGCTa > 2:279344/1‑90 (MQ=255) aCCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCGTTGTTGTCGCTGGGGCTATTACGGCGCGGCt > 2:478848/1‑90 (MQ=255) gcagcagcaAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTaa < 1:126995/90‑1 (MQ=255) ggAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTgcgc < 1:699523/90‑1 (MQ=255) ggAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTgcgc < 1:498016/90‑1 (MQ=255) tATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGcccc < 2:424257/90‑1 (MQ=255) tggtATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTggg < 2:694218/90‑1 (MQ=255) | GGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGG > NZ_CP009273/2435974‑2436147 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |