Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,480,646 | 2 bp→AT | 100% | intergenic (+21/+103) | zwf → / ← gndA | glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,480,646 | 0 | T | A | 100.0% | 118.6 / NA | 36 | intergenic (+21/+104) | zwf/gndA | glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (16/20); total (16/20) | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch. | |||||||||||
* | BS168_reference_genome | 2,480,647 | 0 | A | T | 99.4% | 116.1 / ‑6.8 | 36 | intergenic (+22/+103) | zwf/gndA | glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base T (16/19); total (16/20) | |||||||||||
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 polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
AACAAAACCTTATCTCCTAACTACGAATCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAATAGTTTACTAAGCATTTTCGGTCTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCATCCATTCAGTAT > BS168_reference_genome/2480527‑2480769 || aacaaaACCTTATCACCTAACTACGACTCAGGCACAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAACAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCgg < 7:184519/139‑1 (MQ=255) aacaaaACCTTAGCGCCTACCTCCGAATCCGGCTCAATGGGACCGAAAGAAACTGATGATCTTTTGGTGAAAGACGGCTAACACTGGTGGAACATATAATAAGAAGAAAAACGCCGAAAATGTTTAC‑AAGCATTTTCgg < 7:197786/139‑1 (MQ=255) tATCTCCTAACTACGAATCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACg > 4:209773/1‑139 (MQ=255) aTCTCCTAACTACGAATCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACg < 3:211721/138‑1 (MQ=255) aTCTCCTAACTACGAATCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACg < 3:45726/138‑1 (MQ=255) tAACTACGAATCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGCTTAC‑CAGCATTTTCGG‑CTTTTTTACGCTGAAc > 4:79847/1‑138 (MQ=255) aaCTACGAATCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACaa < 3:209773/139‑1 (MQ=255) aaCTACGAATCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACaa < 4:267412/139‑1 (MQ=255) aaTCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGtat > 6:272854/1‑138 (MQ=255) gCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAAtc > 5:147217/1‑139 (MQ=255) tGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAAt < 3:190524/132‑1 (MQ=255) tGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAAt > 4:190524/1‑132 (MQ=255) aaaGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCTTTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACaaa < 8:254846/139‑1 (MQ=255) aaaGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACaa > 4:168704/1‑138 (MQ=255) tCTGATGATCTTTTTGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAAcccc > 7:285821/1‑138 (MQ=255) ttgttgttttGGTGTATGTCGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATATCATTCTCCTTAATAACATAACCCCCa < 5:221008/130‑1 (MQ=255) ttttGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCtt > 6:40589/1‑136 (MQ=255) ttGGAGAAAGAAGGCTTCCACTCGTGGAACATAAAATACCCAGAAAAAACCCGAAAAAGTTTAC‑AACCATTTTCGG‑CTTTTTTACGCTGAACAAAGGAATAACAAACTCCTTAATAACAAAACCCCCAAAGCTTaac < 6:272199/137‑2 (MQ=255) ggTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATATCACTCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCt > 3:116801/1‑139 (MQ=255) aaaGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCg < 1:83362/138‑1 (MQ=255) aCGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGt < 5:40589/139‑1 (MQ=255) ggCTTACACTGGTGGAAAATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGttt < 8:285821/139‑1 (MQ=255) cTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTc > 6:56657/1‑138 (MQ=255) cTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTGTCTATGGAt > 7:112505/1‑139 (MQ=255) gtggAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTAc > 3:78270/1‑139 (MQ=255) gtggAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTAc > 7:27213/1‑139 (MQ=255) tggAACATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTATCGG‑CTTTTTTACGCTGAACAAAGGTATTACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACt > 6:75016/1‑139 (MQ=255) aCATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCAt < 4:53655/139‑1 (MQ=255) cATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCggg > 4:275718/1‑117 (MQ=255) cATATAATAAGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCggg < 3:275718/117‑1 (MQ=255) taataaGAAGTAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCATCCAtt < 5:213746/139‑1 (MQ=255) taataaGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCATCCAtt < 5:140823/139‑1 (MQ=255) taataaGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCATCCAtt < 4:78270/139‑1 (MQ=255) taataaGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCATCCAtt < 1:186517/139‑1 (MQ=255) taataaGAAGAAAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCATCCAtt < 8:27213/139‑1 (MQ=255) gaagaaAAAAGCCGAAAATGTTTAC‑AAGCATTTTCGG‑CTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCATCCATTCAGTAt > 5:227576/1‑139 (MQ=255) || AACAAAACCTTATCTCCTAACTACGAATCAGGCTCAATGGGACCGAAAGAATCTGATGATCTTTTGGTGAAAGACGGCTTACACTGGTGGAACATATAATAAGAAGAAAAAAGCCGAAATAGTTTACTAAGCATTTTCGGTCTTTTTTACGCTGAACAAAGGTATAACAATCTCCTTAATAACAAAACCCCCAAAGCTTAAGAGCTTCGGGGGTTTCTATTGATTACTTCATCCATTCAGTAT > BS168_reference_genome/2480527‑2480769 |
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: ‑ |