Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AM260480 | 1,842,550 | 0 | C | G | 22.7% | 47.0 / 11.2 | 22 | A724P (GCG→CCG) | h16_B1632 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (9/8); new base G (3/2); total (12/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCGGACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACCGCGTGATCAGCTCCTGCCTGCCGCGGATCTCGCGCCAGGCGTCCGGCTCAAGTTC > AM260480/1842430‑1842679 | cgcgGACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTgg < 1:69641/134‑1 (MQ=255) ggACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGc < 2:837112/125‑1 (MQ=255) ggACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGc > 1:837112/1‑125 (MQ=255) ccGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGcttgg > 2:539391/1‑102 (MQ=255) ccGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGcttgg < 1:539391/106‑5 (MQ=255) ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCgg > 1:212907/1‑133 (MQ=255) ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTgg > 1:775512/1‑104 (MQ=255) ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTgg < 2:775512/104‑1 (MQ=255) ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGc > 1:46967/1‑105 (MQ=255) ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGc < 2:46967/105‑1 (MQ=255) cGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGtt < 1:625724/120‑1 (MQ=255) cGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGtt > 2:625724/1‑120 (MQ=255) ccAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACCGCGTGATCAGCTcctgcct > 1:1018736/1‑134 (MQ=255) cTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGcc < 2:108874/73‑1 (MQ=255) cTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGcc > 1:108874/1‑73 (MQ=255) cTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACcgcg > 1:667084/1‑92 (MQ=255) cTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACcgcg < 2:667084/92‑1 (MQ=255) tctcCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGcc < 2:351669/71‑1 (MQ=255) tctcCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGcc > 1:351669/1‑71 (MQ=255) cGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACcgcg < 2:400049/86‑1 (MQ=255) cGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACcgcg > 1:400049/1‑86 (MQ=255) aGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACCGCGTGATCAGCTCCTGCCTGCCGCGGATCTCGCGCCAGGCGTCCGGCTCAAGTTc > 1:77136/1‑134 (MQ=255) | CGCGGACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACCGCGTGATCAGCTCCTGCCTGCCGCGGATCTCGCGCCAGGCGTCCGGCTCAAGTTC > AM260480/1842430‑1842679 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |