Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,940,921 | 0 | C | . | 64.0% | 19.9 / 34.3 | 25 | noncoding (66/77 nt) | BW25113_RS14640 | tRNA‑Met |
Reads supporting (aligned to +/- strand): ref base C (3/6); new base . (7/9); total (10/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.91e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAA > NZ_CP009273/2940805‑2941051 | aTCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTtgat < 1:498417/139‑1 (MQ=255) aTCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTtgat < 1:62288/139‑1 (MQ=255) gCCCCAGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGCCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAg > 1:501455/1‑138 (MQ=255) gCCCCAGCAACCAAATAAAATTTGATGAAGTAAAGCAGAACGATGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTTAAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGt < 2:501455/139‑1 (MQ=17) aTTTGATGAAGACAGGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg < 2:269321/139‑1 (MQ=255) atgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgggg < 1:333115/139‑1 (MQ=255) atgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTAAAAGc > 1:479096/1‑121 (MQ=255) atgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATGTGATGAAGTAAAAGc < 2:479096/121‑1 (MQ=255) aaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGagca < 1:414819/139‑1 (MQ=255) tACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg > 1:65854/1‑138 (MQ=255) tACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTCAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg > 2:484108/1‑137 (MQ=255) tACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg > 2:418225/1‑137 (MQ=255) tACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTa < 1:484108/139‑1 (MQ=255) tACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTa < 1:418225/139‑1 (MQ=255) cGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat < 2:118021/95‑1 (MQ=255) cGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat > 1:118021/1‑95 (MQ=255) aCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTTGTGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgt < 1:82985/139‑1 (MQ=255) ggTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCTGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGag > 2:93842/1‑117 (MQ=255) ggTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGag < 1:93842/117‑1 (MQ=255) cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGt > 2:416044/1‑68 (MQ=255) cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGt < 1:416044/68‑1 (MQ=255) cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGagca < 2:323213/98‑1 (MQ=255) cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGagca > 1:323213/1‑98 (MQ=255) tcgGTTCAAATCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAat > 2:432684/1‑138 (MQ=255) aaaTCCGG‑CCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGaa < 1:432684/138‑1 (MQ=255) | ATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAA > NZ_CP009273/2940805‑2941051 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |