Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 1,123,238 | 0 | C | A | 59.1% | 6.8 / 18.3 | 24 | intergenic (‑240/+64) | gcvH_1/CJ019_01062 | Glycine cleavage system H protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (2/7); new base A (13/0); total (17/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.11e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑CGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGA‑GGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCA > ACB143/1123117‑1123363 | aaaTTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCGACGACGtgg > 3:22501/1‑132 (MQ=255) aaTTTACGTAGAAGCCGTATTTATTGTGTGGAATGATATCGTAACCGGTGCCTCTGAGAGGAAAAGGGTTTCAGCAACAACCCCTTGTGCCCCGACGGCGAGCAGGCTGTGTAGCGTTTG‑CGCTACGACGTCGc < 2:40789/134‑1 (MQ=255) aaGCCGTATTTATTGTGGGTAATAATAGCGTACCAGGTCCCTCTGAGCGCAAAAGGGTTTGAGCACCAACCCCTTGGGCACCGACGGCGAGGAGGCTGTTTATCGTTTG‑CGCTACCACTTCGCGAAAACGATACa < 8:113181/135‑1 (MQ=255) ccGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCAACCCCCGAGCGAAAACGATACAtgt > 7:218174/1‑135 (MQ=255) acGTATTTATTGTGTGGAATGATATCGTTACAGTTGACTCTGAGAGTAAAATGTTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑ATAAAAGAAGTATCGAAAAAGATACAttt > 6:317039/2‑133 (MQ=255) gtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑CGCTAAGAAGTAGCGCAAACGATACATGTGGTGCCTGGACTg > 7:78019/1‑135 (MQ=255) gtgGAATGATATCGTTACAGGTGCCTCTGAGAGCAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTG‑CGCTACGACGTAGCTAAAACGATACATGTGGTGCCTGGACTg < 8:93592/135‑1 (MQ=255) ggAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCTACGACGTAGCGAAAACGAAAAATGTGGTGCCTGGGCTTGc > 3:302075/1‑135 (MQ=255) gAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCTACGACGTAGCGAAAAAGATAAATTTTGTGCCTGGACTGGCg > 7:70597/1‑135 (MQ=255) aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCTTAGACGTAGAGAAAACGAAACATGT‑GTGCCGTGGCTGGCGA‑g > 7:75905/1‑135 (MQ=255) aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCTAAGACGTAGCGAAAACGATAAAGGTGGTGCCTGGACTGGCGa > 1:120269/1‑135 (MQ=255) tGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCTACGACGTAGCGAAAACGATACATGTGGTGCCAGGACTGGCGC‑gg > 5:76588/1‑135 (MQ=255) tatCGTTACCGGTGCCTCTGCTAGGAAACGGGTTTGCCCATCAACCCCTTGTGCACCGACGGCGCGGCGGCTGTGTATCGTTTT‑CGCTACGACGTCGCTAAAACGATACATGTGGTGCCTGGACTGGCGA‑GGTg < 2:116031/134‑1 (MQ=255) cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑GGCGTCGGCGTAGCGTAAACGCATCATGTGGTGGCTGGACGGGGGA‑GGGGGCaa > 8:316346/1‑135 (MQ=255) cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑CGCTACGCCCTATCGAAAACGATACAAGTGGTTCCTGGACAGGCGA‑GGTGCCAAAGGcc > 1:219023/1‑135 (MQ=255) cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACAGACGGCGAGGAGGCTGTGTATCGTTTT‑AGATTAATCGTTGCGTCAAAGATAACTGTGGTGCCTGGACTTGGGC‑TGTGGCAATGGcc > 1:328055/1‑135 (MQ=255) tgggcccTCTCAGAGCAAACGGGTTTGACCATCAACCCCTGGTGCACCGCCGGCGTGGCGGCTGTGTAGCGTTTG‑CGCTCCCACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGA‑GGTGCCAATGGCCt < 4:302075/130‑1 (MQ=255) gCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑ATCTACGCCCTATCGAAAACCATAAAAGTGGGGGCGGGGCAGGGGC‑GGTGGCAAAGGccgtagg > 8:45635/1‑130 (MQ=255) taaatGGTTTTGCGAATAAAACCTTTGTGCACCTACGGCGGGGTGGGTTTGT‑TGGTTTTCCGCAACCACTTAGCAAAAACGATACATGTGGTGCCTGGACTGGCGA‑GGTGCCAATGGCCTTCGGTCATCTCAc < 7:47279/132‑1 (MQ=255) aaaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGCCTGGCGA‑GGTGCCAATGGCCGTCGGTCATCTCACtt > 3:252300/1‑134 (MQ=255) aaCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑AGCTACGACGTAGCGAAAAAGCTACATGTGGTGCCGGGACGGGCGA‑GGTGCCAATGGCCGTAGGTCATCTACATTTGAGCAGATGATACAGcc > 2:219193/1‑133 (MQ=255) gaggagGCTGTGTATCGTTTT‑AGCGAAGACCGAGCGCAAAAGATAAAAGTGGTGGCGGGACTGGCGAGGGTGCAAAT‑GCCGTCGGTGCTCGGACATTAAGCAGATGATCCACAGAAAACCCCCCGGCGccgccgg > 6:82741/1‑130 (MQ=255) cTGTGTATCGTTTT‑CGTTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGA‑GGTGCCAATGGCCTTCGGTCATCTCACTTTCCGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAg < 4:241219/135‑1 (MQ=255) tATCGTTTT‑GGCGTCGGCGTAGCGGAAACGGTTCATGTGGTGGCGGGGCTGGCGC‑GGTGGCAATGGGCGTCGGTGATCTCCCGTTGAGCAGCTGGTCCAGCGGAAACGCGCCGGCGCCTCGGGCAACCCGGggcc > 8:200650/1‑134 (MQ=255) | AAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTT‑CGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGA‑GGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCA > ACB143/1123117‑1123363 |
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: ‑ |