| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB122 | 2,834,922 | 0 | G | T | 32.0% | 11.5 / 11.4 | 25 | intergenic (‑44/+208) | yqhD/CJ019_02688 | Alcohol dehydrogenase YqhD/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (5/12); new base T (8/0); total (13/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.65e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAATCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGA > ACB122/2834800‑2835046 | gAATCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAACTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTATACCAGTCGTATAATTTTTTTTTTAAAtttttt > 2:124401/1‑135 (MQ=255) aGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTTT‑TTTTTTTTGTGTAATTTTTAttt > 1:199073/1‑132 (MQ=255) gCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAg < 5:183850/135‑1 (MQ=255) aGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCTTAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑TTTTTTTTGGGAAAATTTTATTTAAGAAGCtcgtc > 4:164217/1‑136 (MQ=255) aGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAA‑ATTTTTTTGTGAAATTTTTATTGGAGAAGCgcgtc > 1:139479/1‑136 (MQ=255) cGGGCACTGGTTTGGCTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGgc < 1:130847/136‑1 (MQ=255) cGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGgc < 2:96550/136‑1 (MQ=255) gggCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTTGTTAAATTTTTTTTTAAGAATCACGTCGgcc > 2:75929/1‑135 (MQ=255) aTGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGGAGAAGCTCGTCGgcgc‑gcct > 7:10606/1‑134 (MQ=255) cGAAGATAATTCGCGTGGGGTTGTGTAAGTAAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGAACATTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATa < 3:221359/136‑1 (MQ=255) aGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑AATTTTTTGTTAAATTTTTATTTCAGAATCTCGTCGGCGC‑GCCTTGGGCCCCAAATAga > 6:193356/1‑134 (MQ=255) ggggTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGTAATTTTTATTGCAGAAGCGCGTCGGCGC‑GCCTTGCGCCCCAAATAGCTGAA‑GATTGGtttg > 4:13485/1‑134 (MQ=255) gTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGaa < 5:603/135‑1 (MQ=255) gTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGAAGCTTGCTTCt < 1:75929/136‑1 (MQ=255) gcAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGAAGCTTGCTTCt < 5:220690/134‑1 (MQ=255) aaaaTTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGAAGCTTGCTTCTaaa < 1:196458/136‑1 (MQ=255) gCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGAAGCTTGCTTCTAAAGGCAACCa < 4:55325/136‑1 (MQ=255) aCTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTTGTGTAATTTTTATTGGAGAAACTACTACGCGCGGCATT‑CGGGCCCAATAGCTTTA‑GATTGGTTATGAGCTTGGGTCTAAAGGCCACCCAGTTGGGGTgggcg > 8:199091/1‑133 (MQ=255) tCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCa < 5:194454/133‑1 (MQ=255) tCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCAtt < 5:209045/135‑1 (MQ=255) ttAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAAAAGCGGAA‑GATTTGTTTGAAGCTTGCTTCTAAAAGCAAACAAGTTTGTGTGGCATTCGt > 6:99254/1‑135 (MQ=255) aaaaGTAGACCAGTCGTTT‑ATTATTTTTTGTAATTTTTTTTGAAGAAGATCGTCGGCGC‑GCCTTGGGCACCCAAAAGCTGAAGGATTGTTTAATACATTG‑TTCTTAATGCAACCAATTATGTGTGGGATTAGtgg > 6:182917/1‑133 (MQ=255) aGTAGACCAGTCGTTT‑ATTTTTTGGTGAAATTTTTATTTCAGAAGCTCGTCGGCGC‑GCCTTGCGCACCAAATAGCTGAA‑GATTTGTTATAAGCTTGCTTATAAAGGCCACCAACTGGGGGTGGCTTTCGTATGcct > 8:109470/1‑134 (MQ=255) ccAGTCGGAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCTTTGCGCACCATATAGGTTAA‑GATTGGTTAGAAGCTTGCGTCTAAAGGCAACAAAGTTGGTGTGGCATTCGTCTGCATTCCTgg < 8:118565/136‑1 (MQ=255) tat‑aTTTTTTGGTGAAATTTTTTTTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTATAAGCTTGCGTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTTGATTCCTGGTACCGa > 6:231397/1‑135 (MQ=255) | GAATCAGCTCATCCAGCCGGGCAATGGTTTGGTTACCGAAGATAATTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTAT‑ATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGC‑GCATTGCGCACCAAATAGCTGAA‑GATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGA > ACB122/2834800‑2835046 |
| 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: ‑ |