| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,727,003 | 1 | . | C | 44.4% | ‑6.9 / 16.9 | 9 | noncoding (2182/2904 nt) | rrlG | 23S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base . (2/3); new base C (3/1); total (5/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.24e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTA‑AAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCAAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGC > NC_000913/2726918‑2727139 | ttgggaggaGACCGCCCCAGTCAAACTACCCACCAGGCACTGTCCGCAACCCCGATAAGGGGTCAACGTTAGAACATC‑AACACTACAAGGGTGGTATTTCAAGGACGGCTCCACGAATACTGGCGTACTCGTTTCAAAGcc > 1:285495/4‑141 (MQ=1)ttgggaggaGACCGCCCCAGTCAAACTACCCACCAGGCACTGTCCGCAACCCCGATAAGGGGTCAACGTTAGAACATC‑AACACTACAAGGGTGGTATTTCAAGGACGGCTCCACGAATACTGGCGTACTCGTTTCAAAGcc > 1:372530/4‑141 (MQ=1)ttgggaggaGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTA‑AAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCAAAGcc > 1:1676223/4‑141 (MQ=1) cccacacAAACTACCCACCACACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTA‑AAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCAAAGCCTCCCACCTATCCTac < 1:1655041/135‑1 (MQ=0) acgaCCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTA‑AAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTACAAAGCCTCCCACCTATCCTACACATCAAgg < 1:1384143/138‑1 (MQ=1) aggcacTGTCCGCAACCCCGATAAGGGGTCAACGTTAGAACATC‑AACACTACAAGGGTGGTATTTCAAGGACGGCTCCACGAATACTGGCGTACTCGTTTCAAAGCCTCCCACCTATCCTACACATGTAGGGTCAATGTTc < 1:601444/138‑1 (MQ=0) cgggccccAACCCGGAATACGGGTCAACGATACAACACCAAACACTC‑AAGGCTGGAATTTCAACGTCGCCTCCCTGCAGACTGGCGTCCTCAATACAAAGCCACCAACCTATCCTACACATCAAGGCTCAATGTTCAgtgt < 2:47128/134‑1 (MQ=0) gataaggggTCAACGTTAGAACATC‑AACACTACAAGGGTGGTATTTCAAGGACGGCTCCACGAATACTGGCGTACTCGTTTCAAAGCCTCCCACCTATCCTACACATGTAGGGTCAATGTTCAGTGCCAAGCTGTAGTaaa > 1:2429157/7‑141 (MQ=0) cagtA‑AAGCGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACAATTCAAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcgc > 2:1134417/4‑141 (MQ=0) | TTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTA‑AAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCAAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGC > NC_000913/2726918‑2727139 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |