| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,728,709 | 1 | . | C | 26.7% | 22.9 / 14.8 | 15 | coding (663/1455 nt) | xylB | xylulokinase |
| Reads supporting (aligned to +/- strand): ref base . (6/5); new base C (2/2); total (8/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATG > NC_000913/3728603‑3728850 | ccGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAAtt > 1:129580/1‑149 (MQ=255) gACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCACCTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCCCCGATTTCGAAACTTCCGTTAACAACGCACCAGTAAttt > 1:73252/1‑149 (MQ=255) gACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAAttt > 1:70699/1‑149 (MQ=255) gTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCt < 4:19917/149‑1 (MQ=255) cATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAg < 3:31453/96‑1 (MQ=255) cATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAg > 4:31453/1‑96 (MQ=255) aaCCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACga < 2:129580/148‑1 (MQ=255) cATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATACCCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGAt < 2:70699/149‑1 (MQ=255) cATTCCCACAACAAATTCACCATGTTCATTGTAGCCAACGCCTGCGAAAAATTGCACCGTCGCCATACCCCCACTATTTCTCAACTTCAGGTAACAAAGCAACAGTAATTTCGCTGACTTAGTATAATGCGGGCATATTGTCAAGAGAt < 2:73252/149‑1 (MQ=255) ttGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACt > 4:23731/1‑149 (MQ=255) cACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCa > 6:27095/1‑149 (MQ=255) cACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGCGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCa > 2:124265/1‑149 (MQ=255) gCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGCGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGc < 1:124265/149‑1 (MQ=255) cATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATg > 7:69784/1‑145 (MQ=255) cATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATg < 8:69784/145‑1 (MQ=255) | CCGACGTCCCCAGCGATAACATTGCCTGATTAGCATCAACCATTCCCACACCAACTGCACCAGCTGCATTGTCGCCACCGCCTGCGACAACTGGCACCGTCGCCATA‑CCCCACGCTTTCGCAACTTCAGGTAACAAAGCACCAGTAATTTCGCTGCCTTCGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATG > NC_000913/3728603‑3728850 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |