Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 2,094,808 | 0 | A | C | 47.6% | 18.4 / 10.4 | 44 | L186R (CTG→CGG) | CJ019_01995 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (11/11); new base C (12/8); total (23/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.51e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG > ACB143/2094686‑2094936 | ccATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCCGCCCGCCGGGCCg > 1:15433/1‑137 (MQ=255) ccATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAACGCCCGGCCCCCCAGGccc > 7:202930/1‑136 (MQ=255) ccATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCGCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGACGGCCAGGCGCGCCTGGGCg > 4:252080/1‑137 (MQ=255) cATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCgg > 6:112593/1‑137 (MQ=255) cGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCAGGCCGGGcgg > 7:288727/1‑135 (MQ=255) cGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCCGCGCGCccggcccggcgg > 1:179099/1‑135 (MQ=255) gggCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAgg > 4:217499/1‑137 (MQ=255) ggCCTTGTCCCGGCACAAGAGATTCGCAAGACTCCTCAGCGGCTCTTCATACCCCTTCGGCGCCTGTTGCTGGCGCAAAT‑TCCCCACGCCGCCGTCGATGGCCACGCGCGCGTGGCCGGGCTGGCCTTTTTGCAggg < 8:184381/137‑1 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCAGGCCGGGCGGGCCGTTTTTCAGGGcc > 1:55226/1‑137 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCGGCCCGCCGGGCCGGGCGGGGCGGTTTGAAGGGcc > 1:119589/1‑137 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCCGGCGGGGCGTTTTGCAGGGcc > 3:23915/1‑137 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCAGGCCGGGGTGGGCGTTTTGGCGGGcc > 1:51448/1‑137 (MQ=255) ctcAAGATCCGTCCGCGGCTCTCAAAACCCCTTCGCCTCCGGCGGCTTGCGAAAAA‑CCTCCACCCCGCGCGCCACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaaga < 2:78019/137‑1 (MQ=255) ctcAAGAGACGGCAGCGCCGCTTCAACCCCCTGGGGCGCCGGCGGCTTCCGCAAAT‑CCTGCACGCCGCGCGCGACGGCCAGGCGCGCGTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaaga < 6:117994/137‑1 (MQ=255) cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCGTTTTGCAGGGCCGGTTGGCCGAGGTTGTTGAAGTAc > 2:203737/1‑137 (MQ=255) cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGGCTTTTTGCAGGGCGTGTTGTCCGAGGTTGTTgaagaaa > 7:288665/1‑136 (MQ=255) cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCACGCCCGCCTGGCCCGGCTGGCCGTTTTGCAGGGCCAGTTTGCCGAGGGTTTTGAAGAAc > 6:171164/1‑137 (MQ=255) gaccgctcctCTTCAAACCCCTTCGGCTCGGGCTGCTTGCGAAAAT‑ACCGCCCGCGGCGCCCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg < 1:38318/133‑1 (MQ=255) tccgctgctCTTCATACACCTTCGGCGCGGGCGGCTTCCCAAAAT‑CCGCCCCGCCGCGCTCGCCGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgc < 1:10818/134‑1 (MQ=255) ccAGCGGCTCTACATCCCCCTGCGGCGCCGGGGGCTT‑CCCAAATAACGGCACGCCGCGCGCGACGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgc < 7:123835/136‑1 (MQ=255) tgctCTTCATACACCTTGGGCTCGGGCTGCTTGCTAAAAT‑ACTCCACGCCGCGCGCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag < 3:134422/137‑1 (MQ=255) tgctCTTAATACACGTTGGGCTCGGGCGGGTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag < 4:23915/137‑1 (MQ=255) ttCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCg > 6:255972/1‑137 (MQ=255) tgcGGCTCCGGCTGCTTGCGAAAAT‑CCTCCCCGCCGCGCGCGACGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCg < 2:15433/135‑1 (MQ=255) ggCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGGAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAt > 2:25735/1‑137 (MQ=255) cgccGGCGGCTTGCGCAAAT‑ACCGCCCGCGGCGCTCGAGGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt < 3:90862/135‑1 (MQ=255) cTCGGGCTGCTTGCGAAAAT‑ACTCCCCGCCGCGCTCCAAGGCCAGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt < 5:298483/137‑1 (MQ=255) cTCGGGCTGCTGGCCCAAAA‑ACTGCCCGCCGCGCTCGACGGCCTGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt < 1:74746/137‑1 (MQ=255) cTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGGCGTTTTGCACGGCCTGTTGGGCGAGGTTGTTGAAGAACACTATATTCCGCAGCACCAGGATATGGCGGCttt > 6:11159/1‑137 (MQ=255) cTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCGGGCTGGCCATTTTGCAGGGCGTGTTGGCGGAGGTTGTTGAAGACATGTAGATTCCGCCGCACAAGTAAATGGCGGAgtt > 6:304253/1‑137 (MQ=255) cTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCACGCGCGCCAGGCCGGGCTGGCCTTTTTTCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt > 1:31258/1‑137 (MQ=255) cTCCGGCGGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt < 2:251664/137‑1 (MQ=255) gtgctTGCGAAAAT‑ACTGCCCGCCGCGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGa < 5:112593/136‑1 (MQ=255) tgctTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGAc > 4:194033/1‑137 (MQ=255) tgctTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGGTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGAc > 8:289687/1‑137 (MQ=255) tgccaaaaT‑ACAGCACGCCGCGCCCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa < 8:205583/132‑1 (MQ=255) tGCGCAAAT‑ACTGCACGCCGCTCTCGACGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa < 8:288727/137‑1 (MQ=255) ggcgcAAAA‑CCTCCACCCCGCCCTCGCAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa < 6:292262/136‑1 (MQ=255) ggcgaAAAT‑ACGGCCCGCCGCGCTCTCCGGCATGGCGCGCGTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa < 4:273060/136‑1 (MQ=255) gccgcGCTCGACGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAg < 6:120444/137‑1 (MQ=255) gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGCACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGc > 8:254578/1‑137 (MQ=255) gcTCGAAGGCCACGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCGGTTGGGCGAGGTTTTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGCGGCTTTCCTTGACCCAATTCTTGGCCTAGACAGCGGGc > 4:49034/1‑137 (MQ=255) cGACGGCCAGGCGCGCCGGGCCGGGCTGTCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 8:276758/137‑1 (MQ=255) cGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 2:55226/137‑1 (MQ=255) | CCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAAT‑ACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG > ACB143/2094686‑2094936 |
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: ‑ |