Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 3,440,346 | 0 | A | C | 45.7% | 22.7 / 10.9 | 35 | A24A (GCA→GCC) | CJ019_03266 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (8/11); new base C (16/0); total (24/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.92e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGC > ACB143/3440220‑3440471 | gCAACGGTCTGGCGGCTCGCCAGGCCCTCTACACGCCTGCGCTCGCGGTGCTGCTGTGCGCGCGGGGGCCCCCCGCCGGGGCGCCGGGGCGCGGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGgtg < 8:220117/136‑1 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGgtg > 3:368715/1‑136 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGgtg > 8:108816/1‑136 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGgtg > 4:237919/1‑136 (MQ=255) tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCAgg > 8:101852/1‑136 (MQ=255) cgggaggcGCGCAAGGCCCTCTACACGCGCGCGCGCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCGGGGGCGCTGCTGCTGGCCGGCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc < 1:109158/134‑1 (MQ=255) gggaggaGCGCAAGGCCCTCTACCCGCGCGCGCGCGTGGTGCTGCAGGGCGCGCTGGGGCCACCCGCCGGGGCGCCGTGGCGCGGCGGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGcc < 8:43002/135‑1 (MQ=255) aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCAGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCAc > 7:303069/1‑136 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCTGGCGGCGCCCCGGCCGCCTCGGGGTTGGAGGGCACCGccct > 2:365452/1‑132 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGCCACCGCCGATaa > 1:271519/1‑136 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCACCCCTATaa > 1:273858/1‑136 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCCGGCGGCGCCCCGGCCGCCTAGGTGTGGCAGGCCCACCCCCATaa > 6:15393/1‑136 (MQ=255) gcgcgcgcTCGAGGGGCTGCAGGGCGCGCGTGGGCCACCCGCCGGGGCGCCGGGGCGCTGCTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTg < 8:48649/136‑1 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCCCCGGGGTGTGGCCGGCCCCCGCCGATAATGCGCCGCTCGCTGTTAAgg > 2:43999/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCCCCCCGGCCGCCGGGGTGTGGCAGGCCACCGCCGATAAGTCGCAGCGCGCGGATCAgg > 1:164783/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGGGTGGCCGGCCCCCGCCTATAAGGCGCAGCGCGCGGATCAGGcc > 3:298449/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGGAGGCCAACGCCGAAAATGCGCAGCTCGCGGCTCAGGcc > 3:283915/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGGAGGCCACCGCCGATAAGTCGCCGCTCGCTGATCAGGcc > 7:130206/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCGATAATGCGCAGCTCGCGGATCAGGcc > 3:227235/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCTGGCGGCGCCCGGGCCCCCGGGGGGTGGGAGGCCCCCGCCTATAAAGCGCCGCTCGCGGATCAGGcc > 3:262133/1‑136 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCCCGCCTATAAGGCGCCGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt > 7:291367/1‑136 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCCCGCCGATAAGGCGCCGCTCGCGGCTCAGGCCGATGACGTCGGCCAGCAGt > 3:31444/1‑136 (MQ=255) ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCCGGCGCCCGCCGATAAGGCGCAGCGCGCGGATCTGGCCGAGGCCTTCGGGAGGCAGTTg > 5:26988/1‑136 (MQ=255) ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGGGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg > 7:189168/1‑134 (MQ=255) ggcccGTGGCGGGGGTGCTCGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 8:265412/133‑1 (MQ=255) gcgcgGTGGCGCTGGTGGTCGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 4:177707/135‑1 (MQ=255) cGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgcg < 7:16844/136‑1 (MQ=255) ggCGCTGCTGCTGGCCGACAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGt > 7:89984/1‑136 (MQ=255) ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa < 6:121672/136‑1 (MQ=255) cgccCGCGCTGGCGGCGCCCGGGCACCCGGGGTGTGGCAGGCCAACGCCGATAATGCGCAGCTCGCTGATCAGGCCGATGACGTAGGGAAGCCGTTGGCGGCGAAGGATAGCGCGCACCGGCGGGAGCGGGAGAAg > 2:84876/1‑136 (MQ=255) gccACCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCTATAATGCGCAGCTCGCGGATCAGGCCGCTGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCCCGGTCTGGTGCGGGAGAAgg > 4:160643/1‑136 (MQ=255) ccgccgcCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagc < 8:316119/132‑1 (MQ=255) cggcgCCCGGGCAGCCTGTGGGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTTGAGCGGGAGAAGGCTgcagcagca < 6:111578/136‑1 (MQ=255) cccgcCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGcagcagcag < 2:164783/134‑1 (MQ=255) gcgcCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCagcagcagc > 5:185778/1‑136 (MQ=255) | GCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGC > ACB143/3440220‑3440471 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |