| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 838,331 | 0 | A | C | 53.5% | 25.2 / 43.5 | 43 | T182P (ACC→CCC) | rlmF | 23S rRNA (adenine(1618)‑N(6))‑methyltransferase RlmF |
| Reads supporting (aligned to +/- strand): ref base A (8/12); new base C (0/23); total (8/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.69e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAAC > NZ_CP009273/838242‑838414 | ggCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCAATAAAAACGAGCAATACGACGCGa < 2:1704524/90‑1 (MQ=255) gCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGAc < 2:489448/90‑1 (MQ=255) ttAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCt < 2:662808/90‑1 (MQ=255) ttAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCt < 1:780708/90‑1 (MQ=255) tAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCtt > 1:1617571/1‑90 (MQ=255) gcgcCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTa > 1:803723/1‑90 (MQ=255) aTTCGTTTGCGTCGGCAAAAAGAGAGGGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGCCCTTGTGTAACCCg < 1:483671/90‑1 (MQ=255) ttCGTCTGCGTCGGCAAAAAGAGAGGGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGCCCTTGTGTAACCCGc < 1:1278321/90‑1 (MQ=255) ttCGTCTGCGTCGGCAAAAAGAGAGGGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGc < 1:1479345/90‑1 (MQ=255) tCGTCTGCGTCGGAAAAAAGAGAGGGGGGCGTTTTTTAACGGCATCACCCAAAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGcc < 1:1492329/90‑1 (MQ=255) tCGTCTGCGCGGGCAAAAAGAGAGTGGGGCGTTTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGCCCTTGTGTAACCCGcc < 1:737733/90‑1 (MQ=255) gtttGCGTCGGAAAAAAGAGGGGGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCAt < 1:798461/87‑1 (MQ=255) gcgcCGGAAAAAAGAGAGTGGGGCGATTTTTACCGGCATCTTCCATAAAACCGAGCAATACGACGCGCCTTTGTGTACCCCGCCATTCCa < 1:1850255/90‑1 (MQ=255) gCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCa > 2:1168138/1‑90 (MQ=255) gTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGCCCTTGTGTAACCCGCCATTCCACg < 1:400039/90‑1 (MQ=255) gTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACg > 1:247555/1‑90 (MQ=255) aaaaagGAGAGTGGGGCTTTTTTAACCGGCATCATCCATAAAACCGGGCAATGGGCCGCGACTTTGTGTACCCCGCCATTCCACGATTcc < 1:844043/89‑1 (MQ=255) aaaGAGGGGGGGGCGTTTTTTACCGGCATCTTCAAAAAAAACGACAAATCCGACGCCCCCTTGTTTACCCCGCCATTCCACGATTccgcc < 1:1321667/90‑1 (MQ=255) gggTGATTTTTAACGGCATCAGCCAAAAAAACGAGCAATACGTCGGGCCGTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCAc < 2:973698/90‑1 (MQ=255) gggCGTTTTTTACCGGCATCACCCAAAAAAACGAGCAATACGCCGCGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCAc < 1:453112/90‑1 (MQ=255) gggCGTTTTTTACCGGCATCACCCAAAAAAACGAGAAATACCCGGCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCAc < 1:1386652/90‑1 (MQ=255) gggCGTTTTTTAACGCCACCACCCAAAAAAACGAGCAATGCGCCGCGCCTTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCAc < 1:1572968/90‑1 (MQ=255) gggCGTTTTTAAACGGCATCACCCAAAAAAACGAGAAATCCGCCGCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCAc < 1:1719706/90‑1 (MQ=255) gggCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCAc > 2:730120/1‑90 (MQ=255) gggCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCAc > 2:699229/1‑90 (MQ=255) gggCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCAc < 2:1291105/90‑1 (MQ=255) ttttttaAACGGCATCACCCAAAAAAACGAGAAATACGACGCGCCTTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGc < 1:939173/88‑1 (MQ=255) gttttttACCGGCATCACCCAAAAAACCGAGCAATCCGCCGCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGc < 1:1290133/88‑1 (MQ=255) gttttttAACGGCATCATCCAAAAAAACGAGCAATACGCCGCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGc < 1:699229/88‑1 (MQ=255) cGGCATCATCCAAAAAAACGAGCAATACGACGCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGa < 2:860455/90‑1 (MQ=255) atcCATAAAACCGCGCAATCCGCCGCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTaaa < 1:941028/90‑1 (MQ=255) aaaacCGAGCAATCCGACGCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGt < 1:732122/90‑1 (MQ=255) aaaaaCGAGCAATACGACGCGACCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGt < 2:1693353/90‑1 (MQ=255) aaTCCGCCGCGCCTTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGc < 1:1888441/90‑1 (MQ=255) aaTCCGCCGCGCCCTTGGTTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGc < 1:1273815/90‑1 (MQ=255) atgcgacgCGACCTTGTATACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCt < 1:1508888/87‑1 (MQ=255) atccgacgCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCt < 1:898454/87‑1 (MQ=255) aTACGCCGCGCCCTTGTGTACCCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCt < 1:251303/90‑1 (MQ=255) aTACGACGCGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCt < 2:1090692/90‑1 (MQ=255) aTACGACGCGCCCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCt < 2:153400/90‑1 (MQ=255) aTACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCt < 2:792368/90‑1 (MQ=255) gacgCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAAc > 1:1587818/1‑90 (MQ=255) acgCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGt < 1:739594/74‑1 (MQ=255) acgCGACCTTGTGTAACCCGCCATTCCACGACTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGt > 2:739594/1‑74 (MQ=255) | GGCTTAACCGCGCCATTCGTCTGCGTCGGCAAAAAGAGAGTGGGGCGATTTTTAACGGCATCATCCATAAAAACGAGCAATACGACGCGACCTTGTGTAACCCGCCATTCCACGATTCCGCCGCTGCGGCACGGGCAGGTAGTGAGCGTAAACGCCGTAACCTGGGGCTGAAC > NZ_CP009273/838242‑838414 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |