Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 3,923,486 | 0 | A | C | 51.8% | 1.4 / 19.1 | 27 | W423G (TGG→GGG) | CJ019_03707 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (7/6); new base C (14/0); total (21/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.80e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCACCCCCTGCAAATGCCAAACCCCCAGCGGCTGGTCCG > ACB122/3923357‑3923620 | ggTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAAACAACCCCTGCTCCCCCAGCCCccgcc > 5:183351/1‑135 (MQ=255) cGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGCTCCGCCATCCTGGGGCTGCGAAGCCGACCCGGCAATGTAAAACCAACCCCTGCTCCCCCCGCCCCccccgccc > 7:303892/1‑133 (MQ=255) cGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCCACCCCTGCGCCCCCCGCCCCcgccgccc > 4:244626/1‑133 (MQ=255) cGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAAGCTCACACCCACCCCAGCGGCCGCCGCCCCcgccgccc > 3:23925/1‑133 (MQ=255) ggtgggCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAACCCACCCCCAGCTCCCCCCGCCCCCCCCGCACCcgcctca > 7:402500/1‑134 (MQ=255) ggtgggCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCCACCCCCGCGCCCCCCCCCCCCGCCGCCCGcgcccgt > 2:279787/1‑132 (MQ=255) gggCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCAGACATCCTGGGGCTGCGCAGCAGCCCCGGCCATCTCCAAACAACACCAGGGCCCCCCGCCCCCCGCGCCCGCACCCctgcg > 7:57134/1‑134 (MQ=255) ggATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCAGCTCCCCCAGCCACCGCCGCACCCTCCGCTGCTGTTCa > 7:156779/1‑135 (MQ=255) ccGGGAGTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGAGAAGCAGCCCCGGCAATCTCCAACCCACCCCTGCTCCCCCAGCCCCCGCCGCACGCTCCTCTGCTGTTCACCCGc > 7:279284/1‑135 (MQ=255) aCCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCCCCCCCGGCCATATCCCACCAAACCCCGGTCCCCCCGCCCCCCCCGCACCCGCCAACGCTGGGCACCCCCCAGcccccc > 6:21357/1‑133 (MQ=255) aCCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTGCCCCCGCCCCCGCCGCCCGCACCGCAGCTGGTCACCCGCCAGCCCCGa > 1:94056/1‑135 (MQ=255) tGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCCGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCCGCCCCCGCCGCACGCGCCGCTGCTGTTCACCCGCCACCCCCgccca > 2:92843/1‑131 (MQ=255) tGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCACAAACCAACCCCAGCTCCCCCACCCACCGCCGCACGCTCCTCAGCAGTTCGCCCCCCAGCCCCGCCa > 5:453825/1‑134 (MQ=255) tGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCAGGCAATCTCAAACCAACCCCTGCGCCCCCCCCCCCCCGCGCCCGCTACTGAGCTGTTGACCCGCCAGcccccgca > 8:101549/1‑130 (MQ=255) gCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCGCAAACCAACCCCTGCTGCCCCCCCCCCCCCCGCACGCTCCTCGGCGGTTCCCCCGCACGcccccccaaa > 3:354797/1‑129 (MQ=255) ggACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCCACCCCAGCTCCCCCCGCCACCGCCGCCCGCTCCTCTGGGGGTCACCCGCCAGCCCCGCCCAAAACGCCTCCACCgg > 4:117335/1‑135 (MQ=255) cTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCGCCCCCAGCCACCGCCGCCCGCTCCTCTGCTGTTCACCCCCAAGCCCCGCCAATACCGCCTCACCAGGCCGCCCAcgg > 5:378222/1‑133 (MQ=255) cTGCGAAGCAGCCCCGGCAATCTCAAACCCACCCCGGCTCCCCCACCCACCCCCGCCCGCTCCTCTGCGGGTCACCCCCCAGCCCCGCCAATACCGCCGGAACAGGCCGGCCCCATCGTATTTGTCGCccccccc > 2:20391/1‑134 (MQ=255) ccccccGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGtt < 2:94056/133‑1 (MQ=255) ccccGGCAATCTCAAACCAACCCGTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCAGCAAGCGCCGACAATACCGCCTCATCAGGCGGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTcc < 8:303892/135‑1 (MQ=255) ccccGGCAATCTCAAACCAACCCCTGCTCCCCCCGCCCCCGCCGCACGCTCCTCTGCAGTTCACCCGCAAGCCCCGACAATACCGCCTCCACAGGCCGCCCACTGCGTATTTGTCGCACCCCCCGCGCCCGttac > 8:248594/1‑133 (MQ=255) aTCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCCAGCCCCGACAAAACCCCCTCAACAGGCCGCCCCCAGCGTAATTTTAGCCCCACCGGCGCCAGTTccccccccgg > 5:479544/1‑133 (MQ=255) aTCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCACCCCCTg < 6:41642/135‑1 (MQ=255) aaccaaccCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCACCCCCTGCAAAt < 7:475781/134‑1 (MQ=255) ccaaccCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCACCCCCTGCAAATGcc < 4:354797/135‑1 (MQ=255) cTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCACCAGGCCGCCCCCTGCGTATTTGTCGCCCCACCGGCGCCCGTTCCCCCCCCTGCAAATGCCAACcccc > 8:221432/1‑135 (MQ=255) tCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCACCCCCTGCAAATGCCAAACCCCCAg < 8:402500/134‑1 (MQ=255) aCCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCCCCACCGGCGCCAGTTCCCCCCCCTGCAAATGCCAACCCCCCCGCGGCTGGGccg > 8:145716/1‑135 (MQ=255) | GGTACGCAGCTGGGTGGGCGGGATACCGGGACTTACCTGCAATTGATTGTCGGGACCTTTTTCTGATCCGACATCCTGGGGCTGCGAAGCAGCCCCGGCAATCTCAAACCAACCCCTGCTCCCCCAGCCACCGCCGCACGCTCCTCTGCTGTTCACCCGCAAGCCCCGACAATACCGCCTCAACAGGCCGCCCACTGCGTATTTGTCGCACCACCGGCGCCAGTTCCACCCCCTGCAAATGCCAAACCCCCAGCGGCTGGTCCG > ACB122/3923357‑3923620 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |