Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,097,080:1 | +A | 100% | intergenic (‑104/+46) | azoJ ← / ← yocK | FMN‑dependent NADH‑azoreductase/putative general stress protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,097,080 | 1 | . | A | 96.7% | 98.6 / ‑3.5 | 30 | intergenic (‑104/+46) | azoJ/yocK | FMN‑dependent NADH‑azoreductase/putative general stress protein |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base A (13/16); total (13/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01 | |||||||||||
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. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
CGCTTGATTTTACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑C‑AAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATA > BS168_reference_genome/2096946‑2097202 | cGCTTGATTTTACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTAT‑C‑aaaa < 4:122717/139‑1 (MQ=255) ttttACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑C‑aaaaaaaagc < 4:171505/138‑4 (MQ=255) tACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACATACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACATTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAg < 4:4365/138‑1 (MQ=255) aCAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGt < 4:273003/138‑1 (MQ=255) cAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTa < 4:391329/138‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAatat < 3:49044/139‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAatat < 5:113388/139‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAatat < 7:254972/139‑1 (MQ=255) tGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTActtg > 7:3808/1‑139 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGATTATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCa > 7:380216/1‑139 (MQ=255) cGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtttt < 7:204981/139‑1 (MQ=255) aaaCTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACt > 7:50697/1‑139 (MQ=255) tataTTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCt > 3:290083/1‑139 (MQ=255) tacttacAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCa > 4:334586/1‑139 (MQ=255) acttacATATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCa < 2:359778/138‑1 (MQ=255) tATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGAttt > 7:172730/1‑139 (MQ=255) tGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCa > 5:225384/1‑139 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATcc < 5:89013/122‑1 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATcc > 6:89013/1‑122 (MQ=255) ttttttACTTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAATCAAGTACATATAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTCTTCACGATCAAATATTTCATACGTCTGc > 2:24601/1‑138 (MQ=255) ttttttACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGc > 6:64641/1‑138 (MQ=255) aaGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTg < 5:64641/138‑1 (MQ=255) gTCAATCACATTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATc < 6:223125/139‑1 (MQ=255) aTCACTTTATTTTGAATATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAAtt > 3:184574/1‑139 (MQ=255) ttttGAATATGAGATTTTTTTTAA‑AAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTttg < 7:224153/139‑1 (MQ=255) aaTATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGa > 8:171582/1‑139 (MQ=255) aaTATGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGa > 8:286500/1‑139 (MQ=255) atatGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGAt < 4:391136/139‑1 (MQ=255) tatGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 1:269152/139‑1 (MQ=255) tatGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 3:351493/139‑1 (MQ=255) tatGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 8:172730/139‑1 (MQ=255) tatGAGATTTTTTTT‑CAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 8:380216/139‑1 (MQ=255) | CGCTTGATTTTACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTT‑C‑AAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATA > BS168_reference_genome/2096946‑2097202 |
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: ‑ |