Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | U00096 | 2,947,591 | 0 | A | . | 25.9% | 68.4 / 22.1 | 27 | intergenic (+18/‑16) | metW/metV | tRNA‑Met/tRNA‑Met |
Reads supporting (aligned to +/- strand): ref base A (12/8); new base . (6/1); total (18/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGC > U00096/2947494‑2947737 | ggACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCa < 1:1023238/103‑1 (MQ=255) gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa < 1:1023246/106‑1 (MQ=255) gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg < 1:1023244/109‑1 (MQ=255) gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg < 1:1023259/116‑1 (MQ=255) gcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGc < 1:1023262/86‑1 (MQ=255) tCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGcc < 2:1023265/72‑1 (MQ=255) cccccGCAACCAATCAAATTTTATGTAGTAAAAGCAGTACGGTGACGCGGGGTGGagcagc < 2:1023264/61‑1 (MQ=255) cccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt > 1:1023300/1‑149 (MQ=255) cccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt > 1:1023370/1‑149 (MQ=255) gCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGt > 2:1023311/1‑149 (MQ=255) aaCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGttt > 2:1023355/1‑149 (MQ=255) cAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTg > 1:1023328/1‑149 (MQ=255) aaTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTgg > 2:1023339/1‑149 (MQ=255) taaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtct > 2:1023283/3‑149 (MQ=255) taaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCAAACGGGTGTTTTTTTGTTTCTGGtct > 2:1023292/3‑149 (MQ=255) aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAcc < 1:1023269/119‑1 (MQ=19) aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAcc > 2:1023269/2‑120 (MQ=19) aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtctc > 2:1023367/1‑149 (MQ=255) aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTCGtctc > 1:1023345/1‑149 (MQ=255) aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTcc > 1:1023346/1‑149 (MQ=255) ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATa > 2:1023273/1‑149 (MQ=255) tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATa > 2:1023315/1‑149 (MQ=255) tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATa > 1:1023351/1‑149 (MQ=255) gatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGaa > 2:1023267/1‑109 (MQ=14) tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaa > 2:1023325/1‑149 (MQ=255) gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTaa > 2:1023276/1‑114 (MQ=255) gTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc < 1:1023266/149‑1 (MQ=255) | TGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGC > U00096/2947494‑2947737 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |