Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 3,440,326 | 0 | A | C | 53.1% | ‑5.8 / 12.3 | 32 | S18R (AGC→CGC) | CJ019_03266 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (13/2); new base C (14/3); total (27/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA > ACB143/3440205‑3440456 | gTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGcc > 3:264182/1‑136 (MQ=255) aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCggg > 3:402594/1‑135 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCCGCCTGGgtg > 4:373020/1‑136 (MQ=255) cAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTggggg > 5:382571/1‑133 (MQ=255) cTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCCCAGGCGGCCCCCGGGCCGCCGGGGTGTGGCAgg > 6:63883/1‑135 (MQ=255) tGGAGGCGCGCCAGGCCCTCTACAAGCGCGCGCTCGTGGTGCTGCAGTGAGCGCTTGGGCCCCCCGCCGGGGGGCGGGGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc < 8:232275/135‑1 (MQ=255) ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCCGGCGGCGGCCGGGCAGCCGGGGTGTGGCAGGCCa > 4:218284/1‑136 (MQ=255) cGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCAGCCTGGGTGTGGCAGGCCAACGcc > 1:149589/1‑134 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCCACGccga > 5:84729/1‑132 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCGGGCCAACGCCTATaa > 4:400473/1‑136 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCCGGCCACCCCCgata > 6:379018/1‑135 (MQ=255) agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCACCCCCGATAAGGCGCAGCGCGCTGa > 2:236654/1‑134 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCCGCGCGGGGATCAgg > 4:98774/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCTCGCTGAGCAgg > 4:386359/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAg > 1:10433/1‑135 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACCCCGATAATGCGCAGCGCGCTGCTCAGGc > 8:118066/1‑135 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc > 1:69819/1‑135 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCACGCCCACCCCTATAAGGCGCCCCTCGCTGATAAGGcc > 4:192211/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCCCCCCGGCAGCCTGGGTGTTGGAGGCCCACCCCAATAAAGCGCAGCTCGCTGATCAGGc > 5:101500/1‑135 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCCGGCCGCCGGGGGGGGGGAGGCCACCCCCGATAATGCGCAGCGCGCGGGTCAGGc > 2:407137/1‑135 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCACCGCTGGCGGCGCCCCGGCAGCCTGGGGGTGGGCGGCCCCCCCCGATAAGGCGCAGCGCGCTGGTCAGGc > 6:171729/1‑135 (MQ=255) gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCTGGCGGCGCCCGGGCCGCCGGGGTGGGGCAGGCCACCGCCGATAATGCGCCGCTCGCGGATCAGGCCg > 5:232555/1‑135 (MQ=255) cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGCAGCCGGGGGGTGGCAGGCCACCCCCTATAAGGCGCCCCGCGGCGATCACGCCGCGGACttag > 8:206246/1‑133 (MQ=255) gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATTAGGCGCAGCTCGCTGATCAGGCCGATGCCTTCGGCaa > 8:46739/1‑133 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGGGGCGCCCGGGCCGCCGGGGGGTGGCACGGCAACCCCTCAAAGGCGCCGCGCGCGGCTCAGGCCGCTGACGTCGGCGAggcc > 6:286988/1‑132 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGGCGGCCCACCCCTCTAAGGCGCCGCGCGCGGCTCAGGCCGCGGCCGTCGGCAAGCAg > 6:152701/1‑135 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCGGCCTGGGTGTGGCAGGCCAACCCCAAAAAGGCGCAGCTCGCGGAGCAGGCCGATGACGTCGGCCCGCAg > 7:327344/1‑135 (MQ=255) cgggggCCCCCGGCGGGGGCGGCGGGGCGGGGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt < 1:407137/132‑1 (MQ=255) cGCCGTGGCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 8:41625/135‑1 (MQ=255) ccGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCGGGGGGTGGCAGGCCAACGCCGATAAGGCGCAGCGCGCGGCTCAGGCCGCTGGCTTCGGCAAGCAGTTGGGGGTGAAGGAGcgcgcgg > 6:198475/1‑134 (MQ=255) ccgccAGCGCCGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa < 4:300364/136‑1 (MQ=255) cgccCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa < 6:133231/135‑1 (MQ=255) | GTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA > ACB143/3440205‑3440456 |
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: ‑ |