Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 3,440,346 | 0 | A | C | 52.9% | 11.3 / 28.4 | 34 | A24A (GCA→GCC) | CJ019_03266 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (8/8); new base C (18/0); total (26/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.09e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA > ACB143/3440211‑3440480 | aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAg > 8:283833/1‑137 (MQ=255) gtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTggg > 5:28165/1‑137 (MQ=255) tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgggg > 6:124000/1‑136 (MQ=255) tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCACCCggggg > 2:51532/1‑136 (MQ=255) ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGccg > 1:221203/1‑135 (MQ=255) aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAcc > 5:240260/1‑136 (MQ=255) aTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGcca > 6:243026/1‑136 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCGGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGGGGCCCACGccgagaa > 1:296640/1‑132 (MQ=255) aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCCACGCCGATAAGGCGCAGCTCg > 1:128213/1‑137 (MQ=255) aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCTGGCGGCGCCCCGGCACCCGGGGGGTGGCAGGCCCACCCCCAAAAGGCGCAGCTCGcg > 3:47856/1‑136 (MQ=255) aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGGGGGCGGCCACCGCCAATAAGGACCAGCTCGcg > 3:47857/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGGGGCGCCCCGGCCGCCGGGGTGTGGCGGGCCACCCCCTATAATGAGCCGCTCGCGGCTAAGG‑c > 3:48136/1‑137 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGG‑c > 7:293935/1‑137 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGCTCAGG‑c > 2:96065/1‑137 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCCGGCAGCCTGGGTGGGGGAGGCCACCGCCTATAAGGCGCAGATCCCCGATAAGG‑c > 6:115075/1‑137 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGCCAACGCCTATAAGGCGCCGCTCGCGGATCAGG‑c > 4:206946/1‑137 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGG‑CCg > 7:30078/1‑137 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCGGGGGTGTGGGAGGCCCCCGCCGAGAAGGCGCCGCTCGCGGAGCAGG‑ccc > 6:282849/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGGGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCCACGCCTATAAGGCGCAGCGCGCTGATCAGG‑CCg > 4:226642/1‑137 (MQ=255) gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGGGCCCGGGGCGCCGGGGGGTGGGCGGCCCACGCCGATAAGGCGCCGCTCGCGGCTGAGG‑CCgcg > 1:244495/1‑135 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCACCGCCGATAAGGCGCCGCGCGCGGATCAGG‑CCGATGACTTCGGCAAGCAGtt > 5:283085/1‑137 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCCGCTCGCGGATCAGG‑CCGCTGACGTCGGCCAGCAGtt > 8:170963/1‑137 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCGGGGGGTGGGCCGCCCACCCCTATAAGGCGCCGCACGCGGCTCAGGCCCGAT‑CCGTCGGCAAGCAggt > 8:271464/1‑135 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCCCCGGGGTGTGGGAGGCCAACGCCGATAAAGCGCAGCGCGCTGCTGCGG‑CCGATGACTTGGGCCAGCAGtt > 2:95807/1‑137 (MQ=255) ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGGAGGGCCACCCCGATAAGGCGCAGCGCGCGGGTCAGG‑CCGCTGCCTTGGGCCAGCAGGTGgggg > 8:278861/1‑137 (MQ=255) cccGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCGCGCGGCGCAGG‑CCGAGGACGTCGTCCAGCAGGTGGCGGAGaaa > 8:160760/1‑136 (MQ=255) gggggggTGGCGCTGGGGCGGCCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 5:243026/131‑1 (MQ=255) ggggCCGCCTGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCTCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga < 5:207121/135‑1 (MQ=255) cgggCCGCCCGCGCGGGCCGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCACGCGCACGGTCTGGAGCGGga < 8:30078/135‑1 (MQ=255) ggggggcggcGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag < 6:283085/133‑1 (MQ=255) ggcgccGCGCGGGCAGGCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagc < 3:38093/137‑1 (MQ=255) cgcCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcgc < 2:304220/137‑1 (MQ=255) gCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCTATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTTGAGCGGGAGAAGGCTTCAGCAGCAGCGCTCGATCa < 1:67988/137‑1 (MQ=255) gCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCa < 7:170963/137‑1 (MQ=255) | ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGG‑CCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA > ACB143/3440211‑3440480 |
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: ‑ |