| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,728,110 | 2 | . | C | 48.6% | 55.8 / 81.7 | 37 | coding (1262/1455 nt) | xylB | xylulokinase |
| Reads supporting (aligned to +/- strand): ref base . (10/9); new base C (8/10); total (18/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.46e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAT > NC_000913/3727982‑3728230 | cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGt > 1:151781/1‑139 (MQ=255) aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGGTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑cccccccccccgta < 2:1508016/139‑7 (MQ=255) tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACgg > 2:104650/1‑139 (MQ=255) tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCACCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg > 1:661201/1‑139 (MQ=255) cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGatat < 2:493569/138‑1 (MQ=255) cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGatat > 1:493569/1‑138 (MQ=255) tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAg < 1:104650/139‑1 (MQ=255) aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAg < 1:588183/139‑1 (MQ=255) aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGatat < 2:1478835/126‑1 (MQ=255) aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCCTCCGGTAATCGAGCTGCTGACCGCTGatat > 1:1478835/1‑126 (MQ=255) ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAg > 1:905567/1‑106 (MQ=255) ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAg < 2:905567/106‑1 (MQ=255) tGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc < 1:834722/139‑1 (MQ=255) cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCCTACGGTCATTGGGCCGCTGAACGCTGCTATCCGcccccgcctcacccc > 2:917386/1‑126 (MQ=255) ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGt < 2:151781/139‑1 (MQ=255) gCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCCGCCTGCTGACCCCTGATACCCGCCAGCATCCGACCCCAg > 2:123781/1‑139 (MQ=255) tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccccccccg > 1:219426/1‑78 (MQ=255) tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccccccccg < 2:219426/81‑4 (MQ=255) aTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcg < 2:1170270/139‑1 (MQ=255) tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCCTACGGTAACCGGGCTGCTAACCGCTGATCCCCGCCAGCATCTGACGCCAGCACTCCCTAcgc > 2:996458/1‑139 (MQ=255) aTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgc > 1:771031/1‑139 (MQ=255) ctctGGATTCGCCGCGATCTCCCCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCgccccgcccc < 2:1218477/139‑1 (MQ=255) ggATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt < 2:1204533/139‑1 (MQ=255) ggATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCCCACGGTAATCGACCCGCTGACCGCTGGTATCCGCCACCACCTGACCCCAGGACTCACTCCCCGCCCCCCCCCCAAt > 1:59957/1‑139 (MQ=255) ttCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccccccccg > 2:1268796/1‑57 (MQ=255) ttCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccccccccg < 1:1268796/60‑4 (MQ=255) gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCCTCCCGTCATCGCACTGCTGGCCGCTGAGATCCGCCAGCGTCTGACGCCAGTCCCCACTCCCCCCCCCGCCCCCCATCCAcct > 2:669220/1‑137 (MQ=255) gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg < 1:1215715/101‑1 (MQ=255) gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg > 2:1215715/1‑101 (MQ=255) gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGa > 1:515768/1‑72 (MQ=255) gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGa < 2:515768/72‑1 (MQ=255) aTCTGCGCCAGCCTTGCTGCACCCAGTGCTGGCCCCACATC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACtt < 2:661201/139‑1 (MQ=255) tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg < 1:669220/139‑1 (MQ=255) cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccc > 1:319827/1‑40 (MQ=255) cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑ccccc < 2:319827/40‑1 (MQ=255) cgcCCAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAg < 1:917386/139‑1 (MQ=255) cccAGTGCTGGCCCCACATCCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGc < 1:138120/139‑1 (MQ=255) cccAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc < 1:1293148/94‑1 (MQ=255) cccAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc > 2:1293148/1‑94 (MQ=255) cccAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt > 2:980213/1‑139 (MQ=255) tGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg > 2:1090311/1‑127 (MQ=255) tGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg < 1:1090311/127‑1 (MQ=255) gCCCCACATCCCCCCCCCCATACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc < 2:744675/67‑1 (MQ=255) gCCCCACATCCCCCCCCCCATACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc > 1:744675/1‑67 (MQ=255) | CGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAT > NC_000913/3727982‑3728230 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |