| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 126,970 | 0 | G | C | 27.3% | 33.5 / 14.8 | 22 | V426L (GTT→CTT) | aceF | pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2 |
| Reads supporting (aligned to +/- strand): ref base G (8/8); new base C (3/3); total (11/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGAT > NC_000913/126839‑127115 | gtATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATCTTACTCACTTCGACaaa < 1:51888/149‑1 (MQ=255) cccTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAAcc > 6:27491/1‑148 (MQ=255) ccTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGa < 6:8643/149‑1 (MQ=255) tGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCACATGTTACTCACTTCGACAAAACCGata > 8:2428/1‑149 (MQ=255) aTGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCAc < 6:38517/149‑1 (MQ=255) ctgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGa > 5:21549/1‑140 (MQ=255) ctgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGa < 6:21549/140‑1 (MQ=255) gTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTg > 8:93255/1‑148 (MQ=255) aaGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGtt < 5:27491/149‑1 (MQ=255) gTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCACATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCg < 7:2428/149‑1 (MQ=255) aagTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGt < 7:93255/149‑1 (MQ=255) aaCTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTg > 6:126722/1‑149 (MQ=255) aTCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTaaa < 4:92628/105‑1 (MQ=255) aTCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTaaa > 3:92628/1‑105 (MQ=255) aGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACTAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCAcc > 5:21375/1‑149 (MQ=255) aGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATCTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAAc < 6:114493/126‑1 (MQ=255) aGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATCTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAAc > 5:114493/1‑126 (MQ=255) tctGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGt > 8:74560/1‑149 (MQ=255) tGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATCTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTaaa < 5:82526/88‑1 (MQ=255) tGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATCTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTaaa > 6:82526/1‑88 (MQ=255) gTAACTGGGTAATGATCCCGCATCTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGtt > 5:54440/1‑149 (MQ=255) cATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGAt < 5:126722/149‑1 (MQ=255) | GTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAAAGCCGTTGCTGCAGCTCTTGAGCAGAT > NC_000913/126839‑127115 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |