| 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 | 97.5% | 136.3 / ‑3.4 | 40 | intergenic (‑104/+46) | azoJ/yocK | FMN‑dependent NADH‑azoreductase/putative general stress protein |
| Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (20/19); total (21/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-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. | |||||||||||
TACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑AAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATAAGTCTTT > BS168_reference_genome/2096956‑2097209 | tACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAg < 6:56806/138‑1 (MQ=255) aTAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATc < 2:69475/139‑1 (MQ=255) aaaCTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTaa < 1:13569/139‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTGAATCACTTTATTTTGAATCTGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAatat < 8:283404/139‑1 (MQ=255) tAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAg > 8:50537/1‑124 (MQ=255) tAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAg < 7:50537/124‑1 (MQ=255) tAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCgg > 4:286868/1‑139 (MQ=255) tAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCgg > 4:291327/1‑139 (MQ=255) tAGACATTTGTATTCCTCCAGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATCAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCgt > 6:301190/1‑138 (MQ=255) ttGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTa < 7:79910/119‑1 (MQ=255) ttGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTa > 8:79910/1‑119 (MQ=255) tGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTActtg > 4:162833/1‑139 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCa > 4:77652/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtt > 8:31067/1‑139 (MQ=255) cGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtttt < 1:77285/139‑1 (MQ=255) ttATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTAAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAAt < 2:184988/138‑1 (MQ=255) cAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATc > 4:258207/1‑138 (MQ=255) tAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTTGCATTTTTACTTCGAATCTTGTTCACGATCCAAGa < 1:134148/139‑1 (MQ=255) tAGTATGTTACTTTCTTTTTTACGTAAAGTAAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGa < 2:166838/139‑1 (MQ=255) tATGTTACTTTCTTTTTTACGTAAAGTAAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCTAATCTTGTTCACGATCCAAGAttt > 7:195577/1‑139 (MQ=255) gTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCa > 8:280040/1‑138 (MQ=255) ttACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAatat > 1:232300/1‑82 (MQ=255) ttACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAatat < 2:232300/82‑1 (MQ=255) ttACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTc < 5:64925/136‑1 (MQ=255) ttACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTc > 6:64925/1‑136 (MQ=255) ttACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATa < 5:83134/139‑1 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTg > 7:181610/1‑139 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTg > 1:101308/1‑139 (MQ=255) tATTTTGAATATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCtt > 6:80642/1‑139 (MQ=255) aaTATGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGa > 1:207891/1‑139 (MQ=255) atatGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGCGCTTTTTTGATCAATTGTTTCTTTGTTGGa > 6:69323/1‑138 (MQ=255) atatGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCTATTGTTTCTTTGTTGGa > 7:131043/1‑138 (MQ=255) tatGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 3:291327/139‑1 (MQ=255) tatGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 3:258207/139‑1 (MQ=255) tatGAGATTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 3:286868/139‑1 (MQ=255) tatGAGATTTTTTATCAAAAAAAAGCAAGTACATCTAATCTCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATTCGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 5:301190/139‑1 (MQ=255) tatGAGATTTTTATTCAAAAAAAAGCAAGTACATCTAATCTCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 1:64673/139‑1 (MQ=255) tgagaaTTTTTTTTCAAAAAAAAGCAAGTTCATCTAATATCGGTTTAGTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 8:6113/133‑1 (MQ=255) gagaTTTTTTTTCAAAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATAAGt > 7:51457/1‑139 (MQ=255) tttttttcA‑AAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATAAGTCttt > 3:153555/1‑138 (MQ=255) | TACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑AAAAAAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATAAGTCTTT > BS168_reference_genome/2096956‑2097209 |
| 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: ‑ |