Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 2,583,316 | 0 | C | A | 32.5% | 58.7 / 11.5 | 40 | I113I (ATC→ATA) | CJ019_02435 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (20/7); new base A (1/12); total (21/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.07e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGGGCGGCACCGAAGGCAAGGTGA > KT2440/2583185‑2583446 | ggCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAAt < 8:261156/137‑1 (MQ=255) cGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGgc > 4:86003/1‑137 (MQ=255) agacaagTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGgccg < 1:446270/132‑1 (MQ=255) aaGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATAAAAATGgccgcc > 6:306109/1‑137 (MQ=255) gagTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTc > 3:94408/1‑137 (MQ=255) cAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCgaaga > 5:298987/1‑137 (MQ=255) aGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGaagaaa > 7:203848/1‑137 (MQ=255) tctaaGCCCGAAAGAATCTTAAAGACGTAACCGTGCAGCCGCTCTGAACCGAGGAGTACTACAAGGAGCCCGAGGAGCCGACTGGCCGCCTGTTCTTCAAAATAAAAATGGCCGCCTCGGGAAAGTCTAATAAAGAc < 6:298987/134‑1 (MQ=255) cccGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTcc > 4:42630/1‑137 (MQ=255) gaagCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGa > 8:369622/1‑137 (MQ=255) gaagCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGa > 6:385300/1‑137 (MQ=255) acgCTCAAAGACATACCCGTGCAGCCGCGCTACCCGTCGTATTACCACACGGCGACCGCGCAGCCTACTGGCCGCCGGTTCTTCAACAAAAAAATGCCCGCCTCGGGCAAGTCGAACAACGACGGCTCCGTCTGGAc < 6:227400/135‑1 (MQ=255) aagGTGAACTACTTAACCGGGCCGCCGCGCTACACCTAGGAGTACGACAAGGAGCCCGCGGAGCCGACGGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCGCGGGCAAGTCTAATAAAGACGGCTCCGTCTGGAc < 8:327602/137‑1 (MQ=255) aagATGAAATCCGTACCCGTGCAGCCGCTCTGCACCGAGGAGGACGACCAGGAGCCCGCGGAGCCGACTGGCCGCCTTTTCTTCAACATAAAAATGGCCGCCTCGGGAAAGTCGAAGAAAGACGGCTCCGTCTGGAc < 8:236147/137‑1 (MQ=255) aaccTGAAAGACGAAACCGCGCCGCCGCTCTCCACCGAGGAGACCCACACGGCGACCGAGGAGCCGACGGGCCGCCTGTTCTTCAACATAAAAATGCCCGCCTCGGGCAAGTCGAACAACGACGGCTCCGTCTGGAc < 2:127389/134‑1 (MQ=255) gAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCa > 4:367242/1‑137 (MQ=255) caaaGACGTAACCGTGCAGCCGCTCTTCACCGTGGAGAACGACAAGGAGCCCGAGGAGCCGAGTGGCCGCCTGTTCTTAAACATAAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCa < 2:385771/136‑1 (MQ=255) aCGAAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGAGTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGAAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAagcca < 2:241235/137‑1 (MQ=255) cGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGt > 3:26941/1‑137 (MQ=255) cgCCGCTCTTCACCGTGGATTCCTACCAGGAGACCGAGTCGCCGACTGGCCGACTTTTCTTCAAAATAAAAATGGCCGCCGCGGGCAAGTCCAAGACCCACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGAt < 8:203848/136‑1 (MQ=255) cgCCGCTATTCACGGCGGATTACCAAAAGCAGCCCTAGGCGCCGACTGGCCGCCTGTTCTTAAACATAAAACGGGCCGCCGCGGGCACGTCCAACAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGAt < 2:416323/136‑1 (MQ=255) aGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGAt > 3:15535/1‑137 (MQ=255) ccGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCg < 8:279157/81‑1 (MQ=255) ccGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCg > 7:279157/1‑81 (MQ=255) ccGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGc > 3:225519/1‑137 (MQ=255) ccGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGc > 8:144598/1‑137 (MQ=255) ccGCTCTTCACCGACGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTAAACATCAAAATGGCAGCCTCGGGCAAGACGAAGAAAGACGGCTCCGTCTCGACCCGCAAGCCAGCCCTGTTCGATGc > 5:407397/1‑137 (MQ=255) agTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTg > 7:252665/1‑137 (MQ=255) gaaccaCAAGGACCCCGAGGCGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGc < 6:21884/132‑1 (MQ=255) cggcgcccGCGGCGCCGACTGGCCGCTTTTTCTTAAACAGAAAAATGCCCGCCTCGGGAAAGTCCACGAACCACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTCCGATGCCAAGGGCAAGCCATTGCTGAAAGTc < 4:77854/131‑1 (MQ=255) gaCCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCgg > 1:43781/1‑137 (MQ=255) gaCCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTAGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCgg > 1:456256/1‑137 (MQ=255) agtGGCCGCGTGTTATTCAAAATAAAAATGGCCGACTCAGCCAAGTCCAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGggcgg < 5:44711/135‑1 (MQ=255) aCTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGggcgg > 3:138961/1‑137 (MQ=255) ggCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGGGCGGCAc > 2:388422/1‑137 (MQ=255) gccTTTTCTTAAAAATAAAAATGGCCCCCTCGGGAAAGTCTAATAACGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGGGCGGCACCGaa < 4:405421/137‑1 (MQ=255) ttcttcAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGGGCGGCACCGAAGGCaa < 5:385300/137‑1 (MQ=255) ttcttcAACATAAAAATGGCCGCCTCGGTCAAGTCGAATAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGGGCGGCACCGAAGGCaa < 4:158440/137‑1 (MQ=255) cAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCCGCCCTGTTCGATGCCAAGGGCAAGCCATTGCGGAAAGTCCCGGACATCTGGGGCGGCACCGAAGGCAACGTgg > 4:432203/1‑136 (MQ=255) aaaCATAAAAATGGCCGCCGCGGGCAAGTCTAAGAACGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGGGCGGCACCGAAGGCAAGGTGa < 4:26941/136‑1 (MQ=255) | GGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAGGAGACCGAGGAGCCGACTGGCCGCCTGTTCTTCAACATCAAAATGGCCGCCTCGGGCAAGTCGAAGAAAGACGGCTCCGTCTGGACCCGCAAGCCAGCCCTGTTCGATGCCAAGGGCAAGCCATTGCTGAAAGTCCCGGACATCTGGGGCGGCACCGAAGGCAAGGTGA > KT2440/2583185‑2583446 |
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: ‑ |