Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,327,042 | 0 | G | T | 66.7% | 29.5 / 32.2 | 30 | D312Y (GAT→TAT) | atoB | acetyl‑CoA acetyltransferase |
Reads supporting (aligned to +/- strand): ref base G (6/4); new base T (6/14); total (12/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.39e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAACCTGGGCTTTGATT > NC_000913/2326966‑2327111 | cGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTAtct > 1:381642/1‑81 (MQ=255) aTTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCAt > 2:381645/1‑81 (MQ=255) aTGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAg > 2:381644/1‑81 (MQ=255) gggTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGc > 2:381648/1‑81 (MQ=255) ggTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCt < 2:381609/81‑1 (MQ=255) tGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAAt < 1:381610/80‑1 (MQ=255) tGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAAt > 1:381643/1‑80 (MQ=255) ggCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAg < 2:381611/81‑1 (MQ=255) gCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGc < 2:381612/81‑1 (MQ=255) cAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCAt < 1:381613/81‑1 (MQ=255) cAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCAt > 2:381646/1‑81 (MQ=255) gTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCAttt < 1:381614/81‑1 (MQ=255) tACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTtg > 1:381649/1‑81 (MQ=255) cTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTtgctg < 1:381615/81‑1 (MQ=255) aaaaGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCt < 1:381616/81‑1 (MQ=255) aaaGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCtt < 1:381618/81‑1 (MQ=255) aaaGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCtt < 2:381617/81‑1 (MQ=255) cGTTACAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCg < 1:381619/81‑1 (MQ=255) ttACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGtt > 2:381647/1‑81 (MQ=255) cAACTGGCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTggg < 2:381620/81‑1 (MQ=255) aaCTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGa < 2:381621/81‑1 (MQ=255) ggCGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGaaaaa > 2:381650/1‑81 (MQ=255) cGGGGCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAAcc < 1:381622/81‑1 (MQ=255) gggCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTTCTGCACAGTACCTTGCCGTTGGGAAAAACCTg > 1:381651/1‑81 (MQ=39) ggCTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAACCTgg < 1:381623/81‑1 (MQ=255) cTGCAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAACCTGGGc < 1:381624/81‑1 (MQ=255) tGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAACCTGGGCt > 2:381653/1‑81 (MQ=255) tGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAACCTGGGCt > 2:381656/1‑81 (MQ=255) cAACTGGCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAACCTGGGCtt < 2:381625/80‑1 (MQ=255) gCGGATATTTATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAACCTGGGCTTTGAtt < 2:381626/79‑1 (MQ=255) | CGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGCATTTGCTGCACAGTTCCTTGCCGTTGGGAAAAACCTGGGCTTTGATT > NC_000913/2326966‑2327111 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |