Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 2,113,895 | 0 | T | G | 35.1% | 47.6 / 10.8 | 39 | R663R (CGT→CGG) | CJ019_02007 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (1/23); new base G (13/0); total (15/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.93e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCC > ACB143/2113779‑2114013 | gtcgCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTa < 7:245890/135‑1 (MQ=255) cTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCa < 7:58666/136‑1 (MQ=255) cTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCa < 6:88841/136‑1 (MQ=255) tGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCGCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCt < 1:52940/136‑1 (MQ=255) aGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGaa < 3:289230/135‑1 (MQ=255) gTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTCAATAGATTTTCTGGTACGCCCGCGGCGGCCAGAAATGGt > 6:132167/1‑136 (MQ=255) taCCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACg < 2:34828/134‑1 (MQ=255) ccGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGGGAATATATTTTGTGGTTCGCCCGCGGCGGCCGGGACTGGTGGTCGGCGGCa > 5:119576/1‑135 (MQ=255) gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGGATATAATTGGTTGTTCGGCCGCGGCGGGCGGGAATGGTTTTCGACCTAAAGCAGtt > 1:356735/1‑136 (MQ=255) gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTATTTTCGGGTACGCCACCCGGGGCCGGGACTGGTGGTAGACGTCAACCAGt > 4:128634/1‑135 (MQ=255) gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATATATGTGGTGGTTCGCCCGCGGCGGGCGGGACTGGTGATGGCCCGCAAGCAGtt > 5:76061/1‑136 (MQ=255) gcagTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGGGAATATATGTGGTGGGTCGCCACCGGCGGCCGGAACTGGGGGTAGACGTGAAGCAGt > 3:88763/1‑135 (MQ=255) agTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTATGTGGTGGTAAGCCAGCCGCGGCCGGAACTGGTGGTCGACGGCAAGCAGTTGc > 7:262138/1‑136 (MQ=255) agTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATATATTTGCTGGTTCGCCCGCGGCGGCCGGGAAAGGTGATCGACGGCAAGCAGTTgg > 1:285099/1‑135 (MQ=255) aaGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTAAATATATGTGGGGGGAACCCCGCGGCGGGCTGAACTGGTGATAGAACTCACACAGtttgccg > 7:97320/1‑131 (MQ=255) gATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGGGATCGACGTCAAGCAGTTGCCCGATGGc < 1:150594/135‑1 (MQ=255) tGCCAACGGGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGCGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGAc < 5:72880/136‑1 (MQ=255) gCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGCAACTATGTGGTGGTTAGCCCGCGGCGGGCGGCACAGGGGTTAGGCGGAAAGCAGTTGGCCCAAGGCGccc > 5:342300/1‑136 (MQ=255) gCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGGTGAATCTATGTTCTGGTTCGCCAGCGGCGGGCCGAACGGGTGGTAGACGGCAAAGAGGTTCCCGATGGCgcc > 3:344481/1‑133 (MQ=255) ccAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCt < 8:30213/136‑1 (MQ=255) aaCGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGCTCAATCTATGTGCTGGTAAGCCAGCGGCGGCCTGGAATGGTGATCGCCGTCAAGCAGTTGCCCGATGGCGACCTg > 2:333434/1‑135 (MQ=255) gTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt < 1:279382/136‑1 (MQ=255) gTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt < 6:22847/136‑1 (MQ=255) gTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt < 6:52453/136‑1 (MQ=255) gTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt < 7:144720/136‑1 (MQ=255) tGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt < 3:180214/135‑1 (MQ=255) aggCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGt < 3:258305/134‑1 (MQ=255) cTACGTGCTCGAGCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACa < 8:279080/136‑1 (MQ=255) cGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACatg < 8:230631/135‑3 (MQ=255) gATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGtt < 1:22211/134‑1 (MQ=255) gATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGtt < 5:186771/134‑1 (MQ=255) aTCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACATTCATTTGACAGGTATGtt < 3:58329/133‑1 (MQ=255) tcaggtcaggGTACAGAACCGGGAGACCAGCGAGTCGCGGTCAATCTATGTGCTGGTTAGCCCGCGGCGGCCTGAACAGGTGATAGACGTCAAGCAGTTGCCCCATGGCGCCCAGCAGTTTACAGGAAAGttact > 4:258305/1‑132 (MQ=255) aggggTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTa < 5:113927/134‑1 (MQ=255) gAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGgccagcca < 6:76061/136‑1 (MQ=255) ggAGACCAGCGAGTCGCGATCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACAcc < 7:291086/136‑1 (MQ=255) gagaCCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACAccc < 5:72568/136‑1 (MQ=255) gagaCCAGCGAGTCGCGGTGAATCTATGTGGTGGTTCGCCCGCGGCGGCCGGGACTGGTGATCGACGTCAAGCACTTTGCCCATTGCGCCCGGCAGTTGACAGGTTTGTTCGATGGGGTAGGGcacccacccacca > 6:113927/1‑135 (MQ=255) gagaCCAGCGAGTCGCGGTGAATATATGTGCTGGTTCGCCAGCGGCGGCCTGAACAGGTGGTCGCCGTCAACCAATTTCCCCATGGCGCCCCGCAGTTGGCAGGGATGTTGGATGGGGTTGGCCAGCCACACcccc > 8:13816/1‑136 (MQ=255) | GTCGCTGAGGGCGATGGCAGGCGTTGCACCTATACCGCAGGCAGTAAGGTGGATGCCAACGTGGCCTACGTGCTCGATCAGGTCAGGGTACAGAACCGGGAGACCAGCGAGTCGCGTTCAATCTATGTGCTGGTACGCCAGCGGCGGCCTGAACTGGTGATCGACGTCAAGCAGTTGCCCGATGGCGACCTGCAGTTGACAGGTATGTTCGATGGGGTAGGCCAGCCAGACACCC > ACB143/2113779‑2114013 |
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: ‑ |