Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 3,586,413 | 0 | A | C | 54.2% | 5.8 / 18.5 | 24 | Q131P (CAG→CCG) | CJ019_03399 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (5/6); new base C (12/1); total (17/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.33e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTT‑CGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGAAGCTGG > ACB122/3586285‑3586547 | cGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCCGCCGCCCGCCCCAGCGGGCATT‑AGGCCCGGAccaa > 8:274733/1‑133 (MQ=255) gCACCCTCAATACATCGCACCCAGGGCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTT‑CGGCCAGGACATCGCCGAGCTGCg < 3:138156/136‑1 (MQ=255) aTACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCCGCGCCGCGCCAGAGCGGGATTT‑AGGCCACGCCATAGCCCAGCTGCGGGTTTa > 2:275626/1‑133 (MQ=255) cACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCCCGGCGCCGCCCCCCAGCGGGCTTT‑CCGGCCCGACCTCCCCCCCGAGCGGGATTACCTGGCCCAc > 4:86163/1‑135 (MQ=255) ccccGGGCCGTAGCGCCCGGGGCCGGGGCGGGGCGCCCGGCCCGGCGGGGCCGCCCGGGCGGGCTTGGCGCGGCGCCGCGCCAGAGCGGGATTT‑CGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACaa < 8:224963/135‑1 (MQ=255) aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGCCCGTTCGGGCATTGCGCGGCGCCGCGCCCGCGGGGGGTTT‑GGGGCCGGCCATGGCCGCGGGGCGGGGTTACCTGGCCGCCaa > 6:80859/1‑136 (MQ=255) gggccggggccgCGCCCCCGCCCCGGCCGGGCCGCCGGTCCGGCCAGGGCGCGGCGCCGCGCCAGCGCGGGATTG‑CGGCCAGGACATCGCCGAGCTGCGTGATTCCCTGGCCGACAACGGCAAACGCCGGGTgcg < 2:268049/135‑1 (MQ=255) ctgggccGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGGTTT‑AGGCCCGGACCTCGCCGCGGTGCGGGATTTCCTGGGCGCCAACGGCAAACGCCGGGTGCgcggcg > 7:36759/1‑136 (MQ=255) gcCAGGCACGTCCGGGCCGCCCGTTCGGGCATTGCGCGGCGCCGCGCCAGCGCGGGATTT‑AGGCCCGGGCATCGCCGCGCGGGGGGTTTTCCTGGCCGCCAACCGCCAACGCCCGGTGCGCGGCGGTCCGGGTgcg > 4:26972/1‑136 (MQ=255) gcCAGGCACGTCCGGGCCGCCCGTTCGGGCATTGCGCGGCGCCCCGCCCGGGCGGGCTTT‑CGGGCCGGCCCACGCCGGGCGGCGGGGTTTCCAGGGCGCCCACGGCAAACCCCGGGTGCGCGGCGGTCCCGGTgc > 2:105214/1‑135 (MQ=255) gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGGCCGCGCGGGCTTT‑AGGGCCGGGCATCGCCGGGCTGCGGGCTTACCTGGGCGACCACGGCAAACGCCGGGGGCGCGGCGTTCCGGGTgcg > 8:161703/1‑136 (MQ=255) gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCCGCGCGGGGTTTAGGGCCCGGCCA‑GGCCGCGCGGCGTGGTTTACAGGCCCCCAACGGCCAACGCCCGGGGCGGGGGGGTCCGGGGgcg > 8:200726/1‑136 (MQ=255) gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGCGCGGGGTTT‑AGGCCAGGCCATAGCCGCGCTGCGTGTTTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGGTCCgggggg > 2:250368/1‑132 (MQ=255) gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCCCGCCCGAGCGGGGTTT‑AAGCCCCGGCCTCGCCCCGCGGGGGGATTAACTGGCCCCCAACGGCCAACGCCCGGTGCGCGGGGTGACGGGGgc > 2:15993/1‑135 (MQ=255) gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCCCCCCACCGCGGGAGTT‑TCGCCCCGCCAACCCCGCGCAGCGGGTTTAACAGGGCGCCAACCGCCAACCCCGGGGGCGCCGCCGTCCGGGGgcg > 4:127719/1‑136 (MQ=255) gcCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCCGCGCCCCGCCAGCGCGGGGTTT‑AGGCCCGGGCCTCGCCCCGCGGCGTGATTTCCAGGGCCACAACCGCAAACCCCGGGGGGGCGGCGGACCCGGGgcg > 8:19485/1‑136 (MQ=255) gCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCCGAGCGGGCTTT‑CGGGCCGGACATCGGCGCGCTGCGTGATTACCAGGCCGACCACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCAc > 7:26968/1‑134 (MQ=255) cgccgcgccgCGCCAGACCGGGATTT‑CGGCCACGACATCCCCGAGCTGCGTGATTAACTGGCCGCCCCCCGCAAACCCCGGGTGCGCGGCGTACCCGGTGCGCCCCATGTGCCAGTGTGGATCCTCGGTTCCacc > 8:187016/1‑133 (MQ=255) cAGAGCGGGATTT‑CGGCCAGGACATCGCCGCGCTGCGTGATTACCTGGCCGACACCGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCGGa > 5:109031/1‑136 (MQ=255) cGGGATTT‑CGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTgg > 6:264794/1‑136 (MQ=255) gggTTTT‑CGGCCAGGACAACGCCGAGCTGCGTGATAACCTGGCCGACAACGGCACACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTgg < 1:275626/135‑1 (MQ=255) t‑CGGCCCGGACCTCGCCGAGCTGCGGGATTCCCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGGGCCGGTGTGGATCCTCGGGTCCAGCCTGTTCGGTGCTGACCTGgcggcg < 7:267543/134‑1 (MQ=255) ccAGGACAGCGCCGAGCGGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGAAGCTg < 5:215016/136‑1 (MQ=255) cAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGAAGCTgg < 7:200726/136‑1 (MQ=255) | CGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTT‑CGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGAAGCTGG > ACB122/3586285‑3586547 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |