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.6% | 193.1 / ‑1.3 | 58 | intergenic (‑104/+46) | azoJ/yocK | FMN‑dependent NADH‑azoreductase/putative general stress protein |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base A (25/30); total (26/32) | |||||||||||
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.02e-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. |
TTACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑A‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATA > BS168_reference_genome/2096955‑2097202 | ttACAAATAAAACTGTAGACATTTGTATTGCGCCCGGTAAACATTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATTATATTTTTTTTCAA‑AAAAAAGCAAg > 2:248065/1‑139 (MQ=255) aCAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGt < 6:186015/138‑1 (MQ=255) aCAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGt < 7:142652/138‑1 (MQ=255) aCAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGt < 4:40028/138‑1 (MQ=255) aCAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTa > 2:316395/1‑139 (MQ=255) cAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTa < 2:290799/138‑1 (MQ=255) aTAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATc < 8:264942/139‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAatat < 1:142175/139‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAatat < 8:230647/139‑1 (MQ=255) cTGTAGACATTTGTAATCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAatat < 8:224774/139‑1 (MQ=255) aTTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTAc > 4:157231/1‑139 (MQ=255) ttGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTActt > 5:81959/1‑139 (MQ=255) tGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑a‑aaaa > 1:319210/1‑108 (MQ=255) tGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑a‑aaaa < 2:319210/108‑1 (MQ=255) tGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTGTATTTTGAATATGAGATTATTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTActtg > 4:109259/1‑139 (MQ=255) tATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGGCAATCCGTTTATTTTGATTATGAGATTTTTTTTCAA‑AAAAAAGCATGTACTGCTAATAGCGGGGTActtgat > 5:152482/1‑137 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATc > 5:181788/1‑115 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATc < 6:181788/115‑1 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTAAATCTAATATCGGTGTACTTGCTTGGCa > 2:214945/1‑139 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTAATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCa > 8:236061/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtt > 8:3620/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtt > 1:298660/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtt > 7:36580/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACAACTAATATCGGTGTACTTGCTTGGCAtt > 1:400707/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGAACATCTAATATCGGTGTACTAGCTTGGCttt > 6:210438/1‑139 (MQ=255) ggTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAatat < 4:185555/116‑1 (MQ=255) ggTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAatat > 3:185555/1‑116 (MQ=255) ggTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtttt > 2:301465/1‑138 (MQ=255) gTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑A‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTAc < 8:123232/139‑1 (MQ=255) ttttATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGaa > 2:2425/1‑139 (MQ=255) tttACTTACAAATAGTATGTTACTTTCTTTTTTTCGTAAAGTCAATCACTTTATTTTGAATATGATATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGtt < 3:349118/139‑1 (MQ=255) tacttacAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCa > 5:270462/1‑139 (MQ=255) tacttacAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑ATAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCa > 4:189942/1‑139 (MQ=255) cAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATcc < 1:248065/139‑1 (MQ=255) aaTAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCaa < 5:335760/139‑1 (MQ=255) tAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGa < 1:2425/139‑1 (MQ=255) tAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGa < 7:127194/139‑1 (MQ=255) tATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGATTCTTGTTCACGATCCAAGAttt > 4:198518/1‑139 (MQ=255) tATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGAttt > 4:6165/1‑139 (MQ=255) tATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGAttt > 4:205570/1‑139 (MQ=255) tATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGAttt > 1:394798/1‑139 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtttt > 3:36311/1‑98 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtttt < 4:36311/98‑1 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTg > 3:67799/1‑139 (MQ=255) tttACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAg < 6:270462/138‑1 (MQ=255) cGTAAAGTCAATCACTTTATTTTGAATCTGAGATTTTTTTTCTA‑AAATAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCtttt < 1:214945/139‑1 (MQ=255) cGTAAACTCAATCACTTTATTATGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCtttt < 2:213558/139‑1 (MQ=255) aGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGAt < 3:271858/139‑1 (MQ=255) cTTTATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGttt < 2:394798/139‑1 (MQ=255) tATTTTGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCtt > 5:410274/1‑139 (MQ=255) ttttGAATATGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTttg < 1:316395/139‑1 (MQ=255) tatGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTGATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 2:400707/139‑1 (MQ=255) tatGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 5:309705/139‑1 (MQ=255) tatGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 4:245099/139‑1 (MQ=255) tatGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 7:184471/139‑1 (MQ=255) tatGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 7:236061/139‑1 (MQ=255) tatGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 8:15024/139‑1 (MQ=255) tatGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 1:301465/139‑1 (MQ=255) tatGAGATTTTTTTTCAA‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 8:36580/139‑1 (MQ=255) atGAGATTTTTTTTCAA‑AAAAAAGCTAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCTCGATCGAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 2:146031/138‑1 (MQ=255) | TTACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑A‑AAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATA > BS168_reference_genome/2096955‑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: ‑ |