Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 2,094,808 | 0 | A | C | 62.8% | 6.8 / 11.8 | 47 | L186R (CTG→CGG) | CJ019_01995 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (12/4); new base C (13/14); total (26/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.15e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG > ACB143/2094688‑2094936 | aTGCGGTCCAGTCCCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCCGCTGCTCTTCATACACCTGCGGCGCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCgg < 8:49390/136‑1 (MQ=255) cGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGGCCGGCGCGCcgggccgggcg > 1:352950/1‑135 (MQ=255) gTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGGCAGGCCCGCCGGGCcgggcggg > 8:43323/1‑136 (MQ=255) aGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGGCCGCCTGGCCGGGCAGGcgtg > 4:261618/1‑133 (MQ=255) gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGGTGGGCGTTTTTCAGGGCCGGtt > 8:325828/1‑136 (MQ=255) cccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTg > 4:245376/1‑136 (MQ=255) ccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTAGGCTACGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCGGGCCGTTTTGGCGGGCCGGGTgg > 4:7350/1‑136 (MQ=255) cGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGGGGGGCGTTTTGGAGGGGCtggtggg > 7:70804/1‑135 (MQ=255) cAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGGGCCGTTTTTCCGGGCCGGTTGGGCGAgg > 6:211894/1‑136 (MQ=255) gTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTTGCCGTTTTTCAGGGCCTGGTGGCCGAGgtggtt > 4:129217/1‑136 (MQ=255) tttCTAAATATACGTCAGCTGCTATACAACCACATTCGGCGCGGGTGGCTGGCTCAAATCCTGCACCCCGCGCTCTAAGGCCCGGCTCGCCTGCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt < 8:195162/135‑1 (MQ=255) actcAAGATAGGTCAGCTGCTATTCAACCCCCTTCGGCTCCGGCGGTTTGCCCAAATCCTGCACGCCGCGCGCCAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTga < 5:161604/134‑1 (MQ=255) ccacgATCCGGCAGCGGCGCTGCCTCCACCTGCGGCTCCGGCTGCGTGGGAAAATACTGCCCGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaaga < 2:279988/132‑1 (MQ=255) cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCCCGCCGGGGCCGGCGGGCCGTTTTGCCGGGCCTGTTGGCCGCGGTGGGTgaagca > 4:132009/1‑134 (MQ=255) aCGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGCCGCGCCTGGCCGCGCTGGCCGCTTTGCACGGCCAGGTTGCCGAGGTTGTTTAAGAACTCgat > 1:212744/1‑136 (MQ=255) tccgcggctCATCATACACCTTCGGCTCCGGCTGCTGGCCAAAATACTCCCCGCCGGGCTCGACGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg < 5:156385/130‑1 (MQ=255) cAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCATTTTGCAGGGCCGGTTGGCCGAGGGTGTTGAAGAACTCTATATg > 8:336396/1‑135 (MQ=255) gctgctCTACATCCACCTTCGGCGCGGGCTGCTTGCGAAAATCCTCCCCGCCGCGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcag < 1:44730/136‑1 (MQ=255) ggctCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagca < 1:308143/135‑1 (MQ=255) ctctTCATACCCCTGCGGCGCCGGCTGCTTGCGAAAATACTGCACGCCGCGCGCGAAGGCCTGGCGCGCGTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag < 6:183592/135‑1 (MQ=255) ccacCTGCGGCTCGGGCTGCTGCCGAAAATCGGGCACGCCGCGCGCGAAGGCCCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGat < 1:212339/135‑1 (MQ=255) cgcccccGGCGGCTTCCGAAAATACTCCACGCCGCGCGCTACGGCCTGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc < 2:130779/129‑1 (MQ=255) cGGCTCCGGCGGCTGGCGCAAATACTGCACGCCGCGCGCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCg < 7:254840/135‑1 (MQ=255) ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGa > 3:187609/1‑136 (MQ=255) ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCCCCAGGATATGGCGGa > 4:298911/1‑136 (MQ=255) ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCGGGCTGGCCGTTTTGCAGGGCCTGGTGGCCGAGGGTGTTGAAGAACTCAATATGCAGCAGCCCCAGGATATGGGGGa > 7:160167/1‑136 (MQ=255) cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCCCCTGGCCCGGCGGGGCTTTTTGCAGGGCCAGTTGGCCGGTGTTTGTGTAAAACACTATTATCAACCGCACCAGGATATGGCGGAta > 6:35452/1‑135 (MQ=255) cccGGCGGCTGGCGAAAAACCTCCCCGCCGCGCTCGATGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt < 8:70804/135‑1 (MQ=255) accGGCTGCTTGCGCAAATACTCCCCGCCGCGCTCTAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt < 4:55853/135‑1 (MQ=255) tgcttgctAAAATACTCCACGCCGCGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGa < 1:84703/136‑1 (MQ=255) gctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCTCGGCTGGCCGGGCTGGGCGTTTTTCAGGGCGTGTTGGCCGAGGGTGTTGAAGAACTCTATTTTCAGCAGCACCAGGATATTGCGGAGTTACTTggc > 7:134712/1‑134 (MQ=255) ccgcAAATCCTGCCCGCCGCGCTCCCAGGCCCGGCGCGCGTGGCCGTGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa < 1:362994/135‑1 (MQ=255) cccccgcGCCCGCCGGCCGGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa < 7:319421/133‑1 (MQ=255) aCGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCGGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa > 8:36435/1‑136 (MQ=255) cgccgcGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa > 8:105328/1‑135 (MQ=255) cgccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCGGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAAAGCTTGGCCTCgg > 2:205689/1‑134 (MQ=255) gccgcGCTCGAAGGCCGGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGCCACAATGCTTGGCCTCGa > 5:241131/1‑134 (MQ=255) gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACa > 5:163959/1‑136 (MQ=255) gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTTCAGGGCCTGGTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAAAGCTTGGGCTCGa > 4:154729/1‑134 (MQ=255) gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCGGGCCGTTTTTCAGGGCCTGTTGGCCGCGGTTGTTGAAGAACTCTATATGCCGCCGCACCAGGATATGGCGGAGTTCCTTGACCCAATGCTTGGCCTCGACa > 4:149039/1‑136 (MQ=255) gccgcGCTCGAAGGCCAGGCGCCCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATTCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACa > 4:182278/1‑136 (MQ=255) gcTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCggg > 7:235568/1‑136 (MQ=255) gccggCCCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 3:10928/133‑1 (MQ=255) gacggCCCGGCGGGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 3:30662/133‑1 (MQ=255) gacggCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 3:129632/133‑1 (MQ=255) gacggCCCGGCGCGCCGGCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 7:49333/133‑1 (MQ=255) aGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 2:28748/134‑1 (MQ=255) | ATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG > ACB143/2094688‑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: ‑ |